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ABSTRACT 

We present results from the largest contiguous narrow-band survey in the near-infrared. We 
have used WIRCam/CFHT and the lowOH2 filter (1.187 ± 0.005 pm) to survey «10deg^ 
of contiguous extragalactic sky in the SA22 field. A total of ^ 6000 candidate emission-line 
galaxies are found. We use deep ugrizJK data to obtain robust photometric redshifts. We 
combine our data with the High-redshift Emission Line Survey (HiZELS), explore spectro¬ 
scopic surveys (VVDS, VIPERS) and obtain our own spectroscopic follow-up with KMOS, 
EMOS and MOSEIRE to derive large samples of high-redshift emission-line selected galax¬ 
ies: 3471 Ha emitters at z = 0.8, 1343 [Oin]H-H/3 emitters at z = 1.4 and 572 [Oil] emitters 
at z = 2.2. We probe co-moving volumes of > 10® Mpc® and find significant over-densities, 
including an 8.5 a (spectroscopically confirmed) over-density of Ha emitters at z = 0.81. 
We derive Ha, [Oiii]h-H/3 and [Oil] luminosity functions at z = 0.8,1.4, 2.2, respectively, 
and present implications for future surveys such as Euclid. Our uniquely large volumes/areas 
allow us to sub-divide the samples in thousands of randomised combinations of areas and 
provide a robust empirical measurement of sample/cosmic variance. We show that surveys 
for star-forming/emission-line galaxies at a depth similar to ours can only overcome cosmic- 
variance (errors < 10%) if they are based on volumes > 5 x 10®Mpc®; errors on L* and 
(j)* due to sample (cosmic) variance on surveys probing ^ 10^ Mpc® and ^ 10® Mpc® are 
typically very high: ~ 300% and ^ 40 — 60%, respectively. 

Key words: galaxies: luminosity function, mass function, galaxies: evolution, galaxies: for¬ 
mation, cosmology: observations, cosmology: large-scale structure of Universe, cosmology: 
early Universe 


1 INTRODUCTION 


Determining and understanding the star formation history of the 
Universe is of fundamental importance to improve our understand¬ 
ing of galaxy formation and evolution. Surveys measuring the star 
formation rate density (psfr) as a function of redshift/cosmic time 
suggest that psFR rises up to z ~ 1 — 2 (e.g. Lilly et al.|1996[ Hop-] 


[kins & Beacom|2006[[Karim et al.|2011[|Sobral et al.|2013a[ l and 
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thus reveal that the “epoch” of galaxy formation occurs at z > 1. 
Most interesting, recent studies are also showing that such be¬ 
haviour happens for star-forming galaxies at all masses (e.g.|Karim| 
let al.|201 lllSobral et al.|20T4) . 

There are several star formation indicators that one can use 
to identify and study star-forming galaxies in cosmological vol¬ 
umes. The most direct tracer of recent star formation is the far 
ultra-violet (UV) light, coming directly from very massive, short¬ 
lived stars. Recombination lines resulting from the strong ionis¬ 
ing radiation coming from such stars are also excellent tracers 
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of recent star formation, in particularly Ha, but forbidden lines 
such as [Oll]A3727 can also be used. Studies have also used the 
[Olll] A5007 and other emission lines for this purpose (e.g. |Ly et al.| 
|2007[[Drake et al.|2013| >, although with many caveats such as a po¬ 
tential large AGN contamination. Alternative methods for selecting 
star-forming galaxies are surveying in the far-infrared (FIR) to de¬ 
tect the black-body re-radiation of dust-absorbed UV light from 
young, massive stars, or surveying in the radio for emission of su¬ 
pernova remnants. There is of course a significant difference be¬ 
tween using a star formation indicator to both identify and study 
star-forming galaxies (which can identify and study star-forming 
selected galaxies down to some limit), and using such indicators to 
just measure star formation rates (SFRs) from samples selected by 
some other means, with more complicated and potentially biased 
selection functions. 

Ha stands out as a sensitive star-formation indicator, very 
well-calibrated and not strongly affected by dust extinction in typ¬ 
ical star-forming galaxies, unlike the UV or bluer emission lines 
(typically Ahq = 1 mag, e.g.|Ga m et al.|2010[|Sobral et al.|2012[ 

|Ibar et al.|2013| [Stott et al.|2013at . It is also rnuch more sensitive 

than even the most sensitive FIR or radio surveys, and, with the 
depths that current instrumentation now allow, Ha surveys are able 
to identify both relatively dust-free and dusty star-forming galaxies, 
and thus are ideal for an approximately complete selection of star¬ 
forming galaxies (see e.g. jOteo et al. 12015) . When combined with 
measurements in the FIR, Ha becomes an even better SFR selec¬ 
tor and indicator (e.g. |Garn et al.|2010||Kennicutt et al.|2009||Ibar| 
jet al.|20r3| l. Another potential advantage of Ha is the possibility 
to perform surveys with the narrow-band technique, taking advan¬ 
tage of wide-field optical and near-infrared cameras. Narrow-band 
surveys can probe large areas and, when combined with another 
narrow-band filter at a close-by wavelength, or a broad-band fil¬ 
ter to estimate the continuum, they allow for an effective way to 
obtain clean, complete samples of emission-line selected galaxies 


(e.g. [Bunker et al.' 19951 Ly et al. 

2007[ Geach et al.|2008[ Tadaki 

et al.[[2011 

[Koyama et al. 

2011 

Ly et al.|2011[[Lee et al. 

2012 

Sobral et a 

l.|2012 2013a 

Koyama et al.|2013 [Stroe et al. 

2014 
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The High-redshift(Z) Emission Line Survey (HiZELS; |Geach| 
jet al.|2008{|Sobral et al.|2009a||201^|2013a^ has exploited various 
narrow-band filters in the 2 , J, H and K bands to undertake deep, 
wide surveys for line emitters, with a particular strong emphasis 
on exploring Ha emitters across redshift (e.gdSobral^£talJ20]^ 


[Geach et al.||2012l jSwinbank et al.|[2012a|b| [Stott et al.||2013a|b 

[Sobral et al.|2014^ . HiZELS results in a homogeneous selection of 
Ha emitters at 2 = 0.40, 2 = 0.84, 2 = 1.47 and 2 = 2.23, over 
a few ~ 1 deg^ areas, detecting around 500-1000 emitters at each 
redshift down to limiting Ha SERs of ~ 3 — 10 Mq yr“^ (c.f. 


So- 


[bral et al.[2013a^ . HiZELS allowed for the first fully self-consistent 
measurement of the evolution of the Ha luminosity function from 
2 = 0 to 2 ~ 2, revealing that psFR rises up to 2 ~ 2, using one 
single star formation indicator ( [Sobral et al.|[201^ . The results 
also show that the Ha star formation history can fully reproduce 
the evolution of the stellar mass density since 2 = 2.23 and that 
the typical/characteristic SFR of star-forming galaxies (SFR* ( 2 )) 
has decreased by a factor ~ 13 since that time ( [Sobral et al.|2014] (. 

Results from narrow-band and slitless grism surveys are gen¬ 
erally in good agreement, but some disagreement by factors of a 


few have been reported iLy et al.|2007 

2011[Leeetal.|2012[ So- 

bral et al.|2012l[Drake et al.|2013l[Col 

belt et al.|2013l[Stroe et al. 


sample variance (cosmic variance) due to the relatively small ar¬ 


eas studied. By probing large enough areas, one can overcome the 
effect of cosmic variance (e.g. [Sobral et al.[[2010) , and such data 
can be used to empirically measure the effects of sample variance, 
instead of having to rely on rather uncertain theoretical/model pre¬ 
dictions. Minimising and understanding these discrepancies is not 
only important for our understanding of galaxy formation and evo¬ 
lution, but also as a forecasting tool for upcoming surveys, such 
as Euclid and WFIRST (e.g. [Geach et al.|2010||Wang et al.|2013| 
[Colbert et al.[[20f3) . In particular, determining the bright end of 
the luminosity function with sufficient accuracy, and minimising 
cosmic variance, can have a significant impact on estimates of the 
number of emitters to be recovered with very wide space surveys 
with e.g. Euclid. 


Narrow-band surveys are also sensitive to many other emis¬ 
sion lines besides Ha that can be very useful to (potentially) ex¬ 
tend measurements to higher redshift, but also to provide altema- 
tive/comparative views. These can be used to search for Lya (e.g. 
[Sobral et al.|2009b][Matthee et al.|2014) , but also to trace the evo¬ 
lution of the luminosity function of e.g. oxygen lines, such as [Oil] 
(e.g. jLy et al.jTOQTl [Sobral et al.|201^ . So far studies have mostly 
just covered relatively small areas (e.g. [Ly et al. [20071 [Drake et al.[ 

[2013[ l, probing down to low luminosities, but at the expense of be¬ 
ing strongly affected by cosmic variance, particularly when based 
on a single field. Significant progress requires probing volumes 
~ lOx larger. 


In order to address the current shortcomings/limitations, we 
have undertaken by far the largest narrow-band survey for high 
redshift emission-line galaxies. Observations were made using the 
LowOH2 narrow-band filter on CFHTAVIRCam over a ~ 10 deg^ 
area in the SA22 field and we also combine our new data with sim¬ 
ilar data in the Cosmic Evolution Survey (COSMOS) and the Ultra 
Deep Survey (UDS) from [Sobral et al.[p013a[ . Our study provides 
a major increase in the sample sizes and statistics, but also probes 
a new part of the parameter space. The surveyed area corresponds 
to an area coverage which is ~ 5 — 10 times larger than the largest 
previous emission-line surveys at 2 ~ 1 (e.g. [Sobral et al.|2 0 09^ 
[Ly et al.[2011[[Sobral et al.[2013a^ . We also undertake a significant 
spectroscopic follow-up to improve the robustness of our results. 
Einally, we present the first [Oil] luminosity function at 2 = 2.2. 


The paper is organised in the following way. §2 outlines the 
details of the observations, and describes the data reduction, photo¬ 
metric calibration, source extraction and survey limits. §3 presents 
the narrow-band selection criteria, the final sample of emitters 
and the photometric and spectroscopic redshift analysis, includ¬ 
ing the follow-up observations with MOSEIRE and FMOS. §3 also 
presents the selection of different line emitters within the full sam¬ 
ple, by using colour-colour selections and photometric redshifts. §4 
presents the methods and corrections applied to the data in order to 
derive luminosity functions. Results are presented in §5: the Ha, 
H/3-(-[OlIl], and [Oil] luminosity functions, their evolution to 2 ~ 
2, implications for Euclid, and a quantification of sample (cosmic) 
variance over the full lOdeg^, which includes the HiZELS data. 
Einally, §6 outlines the conclusions. A Hq = 70kms“^ Mpc“^, 
Om = 0.3 and Ha = 0.7 cosmology is used. We use a Salpeter 
IMF and all magnitudes are in the AB system, unless noted other¬ 
wise. 
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Table 1. Observation log for the narrow-band observations conducted with the lowOH2 filter on CFHTAVIRCam. A total of 80 pointings, numbered from 0 
to 79 (Field ID), were obtained with WIRCam, to cover a total area of approximately 10 deg^. The seeing in all observations was in the range 0.5—0.7" and 
conditions were photometric. 


Field ID 

R.A. 

Dec. 

Int. time 

Dates of observations 

3(t limit 


(J2000) 

(J2000) 

(ks/pixel) 


(AB, 2") 

0-28 

22 18 00 to 22 22 00 

-00 04 00 to-1-00 06 24 

1.0 

20-30 Sept, 1-18 Oct, 6 Dec 2011, 4 Oct - 3 Nov 2012 

22.7 

29-53 

22 14 00 to 22 18 00 

-00 04 00 to-I-00 06 24 

1.0 

4-31 Oct, 1-3 Nov 2012 

22.7 

54-79 

22 10 00 to 22 14 00 

-00 04 00 to-I-00 06 24 

1.0 

20 Sept - 6 Dec 2011, 4 Oct - 3 Nov 2012 

22.7 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

Observations were made with the Canada-France-Hawaii Tele¬ 
scope (CFHT) using the Wide-held InfraRed CAMera (WIRCam; 
|Puget et al.|2004|( during three semesters in queued service mode 
(Program IDs 201 IB/029, 20I2A0I9 and 2012B/0I6f] Data were 
obtained using the LowOH2 narrow-band hlter (central wavelength 
of 1187 mij^and FWHM, AA= 10 nm, also referred in this paper as 
NB./), and cover a contiguous region of ~ 10 deg^ (~ 3 x 3 deg^) 
in the SA22 held (centred at 22:16:00, -1-00:18:00; see Table [TJ. 
The SA22 held was chosen for this very wide narrow-band (NB) 
survey as it presents excellent ancillary data over ~ 9.2 deg^ area 
(see Figure[T]for overlap with other surveys). 

Observations were conducted from 20 September 2011 to 6 
December 2011, and were concluded between 4 October 2012 and 

3 November 2012. The seeing over the entire data-set is excel¬ 
lent, in the range 0.5 — 0.7” and observations were done in clear 
conditions. WIRCam’s standard “paw-print” conhguration of four 
2048 X 2048 cryogenically cooled HgCdTe arrays (0.3045"/pixel), 
offset by ~ 1' from each other, probe a region of about 0.11 deg^ 
(21' X 21' held of view at prime focus, see Fig0 at any given time. 
In order to cover the full S A22 held, we obtained 80 different point¬ 
ings. Each pointing was obtained using a dither pattern which re¬ 
sulted in observations being obtained over 10 different positions 
(individual exposures of 100 s) per pointing to hll in the detector 
gaps and minimise the effect of bad pixels and cosmic rays. The 
hnal exposure time per pixel is 1.0ks. A summary of the observa¬ 
tions is given in Table [T] 



Figure 1. Surveyed area in the SA22 with the LowOH2 narrow-band filter 
(NBj) and comparison with other surveys and with the individual WIRCam 
footprint/field-of-view. For Ho emitters at 2 : = 0.81, the surveyed area 
roughly corresponds to a box with ~ 60 X 95 Mpc (physical). The Figure 
also shows (in grey) all naiTow-band detections matched to J and in green 
the location of all line emitters, irrespectively of redshift (see ij^. The over¬ 
lapping regions with CFHTLS W4 (ugriz), UKIDSS DXS (JK), VVDS 
and VIPERS (Zspect) are also shown. 


2.2 Data reduction 

The narrow-band data were reduced with a dedicated pipeline using 
PYTHON and its PYFITS and IMAGE modules, based on PfHiZELS 
( [Sobral et al.|2009^ |2013a[ >. Briefly, we start by median combin¬ 
ing the dark frames to produce master darks and then use them to 
dark subtract the individual science frames. We obtain first-pass 
flat fields by median combining jittered science frames, and use 
those to flatten the data. We then run SEXTRACTOR ( [Bertin &] 
|Amouts||l996[ ) on the first-pass flattened frames to produce indi¬ 
vidual masks. We use those to mask out all individual sources, 
and, excluding each frame that is being flattened, we produce a 
final flat field for that frame and flatten the frame. We then use 


SCAMlj^ I Bertin||2006| to fit a World Coordinate System (WCS) 
by matching sources detected in individual reduced science images 
to those in 2MASS jSkrutskie et al.|2006| l catalogue. We also use 
SCAMP to correct for distortions across the field of view by fitting 
a third order polynomial. Frames are also normalised individually 
to the same zeropoint by computing the ratio between the expected 
flux/magnitude from 2MASS and that found in our data (see more 
details in Section [23) . For both steps we are able to use on aver¬ 
age ~ 75 stars per individual frame, resulting in astrometric solu¬ 
tions with typical rms of ~ 0.15”. Finally, the individual reduced 
frames from the four detectors are median combined using SWARlj^ 
( |Bertin|[2010| ), to obtain stacked reduced data for the entire field. 
For the broadband J (and K) comparison data, we use UKIDSS- 


For the broad band J (and K) comp a 
DXS-DR1 (Q ( [Lawrence et al.|2007t . 


^ FP7/2007-2013: The research leading to these results has received fund¬ 
ing from the European Community’s Seventh Framework Programme 

(FP7/2007-2013) under grant agreement number RG226604 (OPTICON). ^ http://www.astromatic.net/software/scamp 

^ We fully confirm the central wavelength by exploring all spectroscopic ^ http://www.astromatic.net/software/swarp 

redshifts in our study. ^ http://surveys.roe.ac.uk/wsa/ 
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2.3 Source Extraction and Snrvey Limits 

We obtain the magnitude zeropoint (ZP) by comparing the mag¬ 
nitudes of the sources in the 2MASS catalogue to those in our 
dataset, excluding the faintest (J > 17, low S/N in 2MASS) and 
the brightest (J < 12, saturated in our data) sources. In order to 
simplify the analysis, and once accurate ZPs are determined for 
each stacked image, we scale all images such that ZPs are set to 
25, including the broad-band images (J, UKIDSS). Sources were 
extracted using SEXTRACTOR ( [Bertin & Arnouts||1996^ on both 
narrow-band (NB) and broad-band J (BB) images, using 2” aper¬ 
tures. The 3cr AB-magnitude limit for the survey is 22.7 (5 ct: 22.15; 
c.f. [Matthee et al.||2014[ >, corresponding to an emission line flux 
limit of 8 X 10~^' ergs“^ cm~^. This limit is computed by mea¬ 
suring the average background rms of the narrow-band images in 1 
million empty 2" diameter apertures and for point-like sources (no 
aperture correction is applied as it is very small given the excellent 
seeing, assuming sources are point sources); this is the aperture we 
use throughout the paper for all measurements and corresponds to 
about ~ 16 kpc physical diameter for the 3 lines investigated here. 
We note that because we use random aperture measurements, the 
rms that we measure already accounts for correlations in the noise. 


3 SELECTION OE EMITTERS 

Down to the 3 (t limit we detect 346,244 sources in our narrow- 
band images. Once catalogs with sources in the narrow-band (NB) 
and in the broad-band (BB) are made, they are matched using a sky 
algorithm with a maximum separation of 1”. Narrow-band sources 
with no matching broad-band source are likely to be spurious, but 
they are kept in the catalogue and assigned a J upper limit. 


3.1 Emission line candidates 


2013af . Here we correct for this effect using CFHTL!^ 
2 -band, which is the closest band on the blue side of J. Our colour 
correction (empirical) is given by: 

BBj - NBj = {BBj - NBj)o -b 0.04(z' - BBj) - 0.05. (4) 
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For sources with no z band (< 3 a) available (2%, either 
because they are too faint in 2 or because they are masked in 
CFHTLS), we apply the average correction obtained for all the 
sources which have reliable 2 detections (—0.02). Thus, we note 
that this is, on average, a very small correction. 

We classify as potential emitters the sources that have E > 3 
(see Figure]^. Table 2 indicates the major emission lines expected. 
The second criterion for an excess source to be an emitter is that 
the emission line must have an observed-frame equivalent width 
(EW, the ratio of the line flux and the continuum flux densities) 
higher than the scatter at bright magnitudes. This step avoids se¬ 
lecting sources with highly non-uniform continua (with e.g. strong 
features). We compute EWs by using: 


EW = AAnBj -j- 

/bBj 


/nBj - /bBj _ 

/nBj ( AAnBj / AAbbj ) 


(5) 


AAnbj and AAbb j are the widths of the filters and /nbj and /bbj 
are the flux densities for the narrow and broad band respectively. In 
order to identify a source as a potential line emitter we require it 
to have EW (observed) higher than 30 A, corresponding to an ex¬ 
cess of J — NBj > 0.3. Note that this will correspond to different 
rest-frame equivalent widths depending on the line/redshift being 
looked at. For z = 0.81, the rest-frame EW limit (Hq:-i-[Nii]) is 
~ 17 A, while for [OlIl]-l-H(?, at 2 ~ 1.4, and [Oil], at 2 = 2.2, 
it is 12.5 and 9 A, respectively. These are all relatively low equiv¬ 
alent widths for sources at 2 ~ 1 — 2 (e.g. |Fumagalli et al.|2012[ 
[Sobral et al.|2014l ), and thus the sample completeness will be high 
(note that we take this into account when obtaining completeness 
corrections). 

Fluxes of emission lines are calculated as follows: 


In order to robustly select sources that show a real colour excess of 
the narrow-band over the broad-band, instead of just random scatter 
or uncertainty in the measurements, two criteria are used. First, the 
parameter E ( [Bunker et al.|1995t is used to quantify the real excess 
compared to an excess due to random scatter. This means that the 
difference between counts in the narrow-band and the broad-band 
must be higher than the total error times E: 

Cnbj — Cbbj > E5. (1) 


= AAi' 


Inbj - fsBj 


1 — (AAjvbj/AAbsj) 


( 6 ) 


Throughout these calculations, the conversion from magni¬ 
tudes to flux densities (f^) is: 


/„ = 10^^ X IQ-OW^AB+iSS) ^ 


(V) 


where the flux density is given in Jansky and niAB is the magni¬ 
tude in AB. 


Here cnBj and cbBj are the counts in the narrow- and broad-band 
respectively, while S is the total photometric error, the combination 
of the errors in both bands: 



where Tap is the radius of the apertures in pixels and a the RMS 
per pixel in each band. These two equations can be combined into 
the following equation for E jSobral et al.[2013^ : 

1 _ -|n-0.4(BB-NB,) 

s = , - (3) 

ZP is the zeropoint of the NBj (which is the same as the BBj, 
because they are both scaled to ZP = 25 in our analysis). We note 
that the central wavelength of the narrow-band is not perfectly in 
the centre of the broadband (J), but rather at the blue end of the 
filter, similarly to other narrow-band filters in the J band (see e.g. 


Using our selection criteria, out of the 346,244 NB sources in¬ 
dividually detected, 8,599 emitters were selected as potential line 
emitters. However, many of these are still likely to be artefacts 
and/or sources in very noisy regions. This is because these numbers 
do not include any masking for bright stars and their haloes/nor for 
other artefacts (see e.g. [Geach et al.|2008[[Sobral et al.|2009a^ . We 
therefore clean our list of potential emitters by visually inspecting 
all candidates before flagging them as final emitters and produce 
a final mask. After visually checking all the candidates (obtained 
before masking), 2,284 are marked as spurious or artefacts (fully 
consistent with e.g. [Geach et al.|2008^ [Sobral et al.[[2013^ , with 
the vast majority being either fake sources/artefacts caused by the 
many bright (J < 12) stars in our full coverage, or real sources 
which have their NB flux boosted (relative to J) due to enhanced 

® Canada-France-Hawaii Telescope Legacy Survey 
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Table 2. Redshifts at which the LowOH2 filter (Ac = 1187 nm, AA = 10 
nm), targets strong emission lines and corresponding probed co-moving 
volumes per square degree. There are, of course, many other rarer emis¬ 
sion lines which the filter can obtain - see §3.3. * see |Matthee et al.H2014) 
for full details on the search for Lya at z = 8.8. 


Emission line 

^0 

(nm) 

z 

Volume 

10® (Mpc^deg“^) 

Ha 

656.3 

0.81 ± 0.01 

0.11 

[OIII] 

500.7 

1.37 ±0.01 

0.23 

[OIII] 

495.9 

1.39 ±0.01 

0.23 

H/3 

486.1 

1.44 ± 0.01 

0.24 

[On] 

372.7 

2.18 ±0.02 

0.36 

Lyo* 

121.6 

8.76 ± 0.04 

0.52 


stripes, again caused by bright stars. We note that the vast major¬ 
ity of these sources could be masked in an automated way, but 
due to the very wide area of the survey and the different artefacts 
caused by stars with different bright magnitudes, our method (visu¬ 
ally checking every single potential emitter) assured that no region 
would be overlooked. The remaining sources are excluded because 
they are detected in low S/N or noisy regions. This leads to a sam¬ 
ple of 6,315 potential emitters (see Figure[^, covering an effective 
area (after masking) of 7.6 deg^. 

Following |Sobral et al.| ( |201^ , we then exclude potential stars 
within the sample by using an optical vs a near-IR colour. In our 
case, we use g — z versus J — K and identify blue sources in 
J — K which are very red in p — ^ This results in identifying 
339 potential stars, which we remove from the sample, resulting 
in a final sample of robust line emitters of 5976 line emitters. This 
density of emitters is in very good agreement with similar surveys 
of smaller areas ( [Sobral et al.|2013^ . 

We present the catalogue of emitters containing the final 6315 
excess sources, and identify candidate stars with the flag -1 in Ap¬ 
pendix]^ The catalogue contains photometry in narrow-hand J 
(from this work, 2") and J (from UKIDSS, obtained by us, also 
with 2” apertures), line flux. Equivalent Widths and E. Table 3 
summarises the catalogue which is available on-line together with 
the paper. 

The data provides the opportunity to identify large statistical 
samples of emission-line galaxies, mainly Ha at z = 0.81, but 
also [OlIl]-l-H/3 at z ~ 1.3 — 1.4, [Oil] at z = 2.18, and Lyman- 
Q candidates at z = 8.8 (see |Matthee et al.|2014| ). Redshifts and 
probed volumes for the strongest/main emission lines to which our 
survey is sensitive are found in Table]^ 


3.2 Photometry and photometric redshifts 

We combine deep (~ 23 AB) J and K data (UKIDSS DXS DRIO) 
from UKIRTAVFCAM ( [Lawrence et al.|2007) with CFHT Legacy 
survey data in ugriz (limit ~ 26 AB), to produce a photometric 
catalogue down to the 5 ct depth of the J band data (J-selected) us¬ 
ing 2” apertures. We then match the catalogue with our catalogue 


^ We identify stars as sources with (g — z) > 7 X (J — K) + 2 A {g — z > 
1.4) A {J — K) > —0.5. We check that this procedure eliminates the bulk 
of spectroscopic stars in both SA22 and COSMOS NBj samples without 
excluding real line emitters. 



-2 1 . 

14 15 16 17 18 19 20 21 22 23 

NBj (mag, AB) 


Figure 2. Colour-magnitude diagram for our NBj detected sources. The 
dotted horizontal line is for an observed EW of 30 A, which corresponds 
to J-NB J = 0.3. The 3 and 5-S selection curves are shown for the average 
depth of the survey for reference (but note that some small variations exist 
across the survey). Candidate line emitters that passed the EW, S and vi¬ 
sual filtering (to exclude spurious sources or artefacts) are shown in colour, 
while all other non-excess sources are shown in grey. Emitters classed as 
Ha (z = 0.81), [OlII] (-I-H/3) at z = 1.37 and [Oil] at z = 2.2, or unclas¬ 
sified, are plotted in different colours as indicated by the key. Note that, as 
expected, higher redshift emitters are always preferentially found at fainter 
and fainter magnitudes. Thus, for example, a survey a magnitude shallower 
would recover little to no [Oil] emitters at z = 2.2. 

of NBj emitter^ By using the J-selected catalogue, we measure 
PSF matched magnitudes in ugrizK (2 " apertures), and compile 
a ugriz JK catalogue. We use the photometric catalogue to distin¬ 
guish between different line emitters by using colour-colour diag¬ 
nostics, but also by deriving and exploring photometric redshifts. 

We compute photometric redshifts by using EAZY jBrammer| 
jet al.|200^ , which contains a wealth of templates with the main 
emission lines included. We use ugriz JK photometry, but add our 
narrow-band photometry as well (we also run EAZY without the 
NB to control for this addition). In total, we obtain photometric red¬ 
shifts for 5953 emitters (> 99 %). The photometric redshift distri¬ 
bution of our sample of emitters is shown in Eigure]^ We find that 
the EAZY photometric redshift distribution has clear peaks at the 
redshifts of our expected strong emission lines (e.g. Ha at z = 0.8, 
H/3/[Oin] at z = 1.4; see Eigure]^. We also find a tentative peak 
around z ~ 0.25, which may be driven by [SlIl]9530.6A emitters 
(e.g. jMilvang-Jensen et al.|2013| l - we confirm some spectroscop- 
ically (see Eigure]^. We note that the peaks corresponding to the 
redshifts of strong emission lines are clearly enhanced by including 
the narrow-hand data. These, probably more accurate, photometric 
redshifts are used additionally to colour-colour selection and spec¬ 
troscopic redshifts to distinguish between various emission lines, as 
outlined in §3.5. Furthermore, in Figure|^we show how our EAZY 
photometric redshifts (with the inclusion of the NB) compare with 
all the spectroscopic redshifts available for the sample, revealing 
good agreement (the spectroscopic redshift distribution of the sam¬ 
ple with spectroscopic redshifts is shown in Eigure]^. 

® For those without a J band > 5 cr detection (5%) we measure their J 
magnitude centred on the NB j position and, if undetected above 2(7 (8 
sources, 0.13%), we assign the 2 a limit. 
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Figure 3. Left: Photometric redshift distributions obtained with EAZY (including the NB data) for our NBj selected line emitters in the SA22 field and 
comparison with the final samples of the different emitters. This shows that the sample of line emitters is dominated by Ha emitters at z = 0.8, followed by 
H/3+[OlII] emitters (z ~ 1.4), while [Oil] emitters at z ~ 2.2 represent a lower fraction. Right:'The distribution of spectroscopic redshifts of the full sample 
of line emitters, including our follow-up observations with FMOS on Subaru and MOSFIRE on Keck. This clearly shows that the sample of spectroscopically 
confirmed line emitters is dominated by Ha emitters, followed by [OlII], [On], H/3 and [SlI] emitters. Other much rarer lines are found, including Pa7l0870, 
Nll0406,CI9853/9827, [Slll]9533, Hell8237, Ol7774 OII7330/7320, [Arlll]7135, [Ol]6363/6300, Hel5876, [Slll]6311-i-Hell6311. Note that the sample with 
available spectroscopic redshifts is still dominated by optical follow-up (and observed-optical selection), and thus is highly biased against the higher redshift 
line emitters. Dashed lines indicate the redshifts of the major emission lines studied in this paper. 


3.3 Spectroscopic redshifts: literature compilation 

Spectroscopic redshifts are available from the VIMOS-VLT Deep 
Survey (VVDsQ which covered 4 deg^ in the SA22 field. Addi¬ 
tional spectra are available from VIPERS jGarilli et al.|20l4) . In 
total, 289 of our candidate line emitters have a spectroscopic red¬ 
shift from either VVDS or VIPERS. However, and particularly for 
redshifts above 1 (but even for z ~ 0.8), it becomes significantly 
more difficult to identify emission lines in the VVDS survey (the 
survey selection is I<22.5 and there are few and weak lines in the 
optical), as the strongest lines are shifted to the (near)-infrared, or 
still at bluer wavelengths (like Lyman-a). Thus, targeted spectro¬ 
scopic follow-up (directly in the observed NIR) is ideal to increase 
the fraction of spectroscopically confirmed line emitters, and allow 
a more representative evaluation of the range of line emitters within 
the sample that is not biased towards rest-frame UV bright galaxies. 


3.4 Follow-up observations with MOSFIRE and FMOS 


In order to significantly complement spectroscopic redshifts avail¬ 
able from the literature, we follow-up a significant fraction of the 
sources. We obtained spectra for some of our sources using KMOS 
I Sharpies et al. J20]3k and these are presented in |Sobral et al.| 
( 2013b^ and Stott et ai.H2013b[ l. We also followed up some of the 
bright line emitters with WHT and NTT; these are presented in |So-| 
Ibral et al.| ( [M5l >. 

Furthermore, we have observed large samples of our emit- 
ters using FMOS jK imura e t al.|2010^ on Subaru and MOSFIRE 
( [McLean et al.|2010[|2012| > on Keck. Both instruments are shown 
to be tremendously efficient for our targets. This is because not only 
do we know where one of the main emission lines should be found 
(within a low OH emission window, thus maximising the S/N), but 


http://cesam.oamp.fr/vvdsproject/index.html; Le Fevre et al. 


20131 


we also know that the vast majority have fluxes high enough to be 
detected in modest exposure times. 


3.4.1 FMOS observations 

FMOS observations were taken on 15 June 2014, under clear con¬ 
ditions and good seeing (0.7-0.9")- We observed sources within the 
SA22 field in two different configurations, centred on 22 19 57.62 
-bOO 19 35.28 (PI) and 22 19 05.06 -l-OO 52 34.19 (P2). We selected 
a total of 128 sources from our candidate line emitters in SA22 for 
PI and another set of 128 sources from our candidate line emitters 
in SA22 for P2. The remaining fibres (72) per configuration tar¬ 
geted fillers, calibration stars or were disabled. We used the J-Long 
{R — 2200) setting (high resolution mode). Individual exposure 
times were 900s in each nodding position and we obtained two of 
each position for a total exposure time of 3.6ks per pointing, cor¬ 
responding to 1.8 ks of on-target exposure time due to cross-beam 
switching. 

Apart from the raw science frames, we obtained (per configu¬ 
ration) dome-flats and Th-Ar spectral calibration arcs which were 
used for the reduction. In order to reduce the data, we used the 
Subaru FMOS reduction pipeline, FMOS Image-Based REduction 
Package (FIBRE-PAC, |Iwamuro et al.|2012[ l. FIBRE-PAC is a com¬ 
bination of IRAF tasks and C programs using the CFITSIO library 
( |Pence|1999| l. Details are found in jlwamuro et al.| ( |2012| l. Briefly, 
data are flat-fielded and bad pixels removed to begin with. After 
that, corrections are applied to fix spatial and spectral distortions 
present and then wavelength calibration is performed. An initial 
background subtraction is then achieved using the ABAB nodding 
pattern of the telescope to perform an A-B sky subtraction. Further 
bad-pixel, detector cross talk, bias difference, distortion, residual 
background and sky corrections are applied. The 2-d spectra are 
combined, which in cross beam switching mode means inverting 
and adding the negative B spectra to the A spectra. The final step 
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Figure 4. Top: Some examples of our MOSFIRE spectra. Top left, one of the confirmed [OlII] emitters with a strong detection of H/3, where the emission 
line being measured and detected by our narrow-band filter is actually [OlII]4959. Top middle', a very luminous broad-line AGN with extremely broad Ha 
line, implying a very massive black hole of ~ 10® Mq. Top right. A high EW Ho emitter, representative of the less massive population of Ha emitters. The 
locations of the strongest OH lines are indicated as shaded lines, but we note that other weaker OH lines are also present. No smoothing has been applied. 
Middle: Representative examples of spectroscopically confirmed Ho emitters in our sample within our FMOS data-set. No smoothing has been used. The 
strongest OH emission lines are represented with shaded regions (but weaker OH lines are not indicated). Bottom: Representative examples of spectroscopically 
confirmed [OlII] emitters in our sample within our FMOS data-set. No smoothing has been used. 


is an initial flux calibration, which we cross check with our NB 
estimated fluxes (see |Stott et al.|2013b| for more details). 

We obtained redshifts by first identifying the emission line 
within the NB j filter profile wavelength range, and then fitting with 
various combinations. For most spectroscopic redshift determina¬ 
tions there were at least 2 lines, but whenever only one line was 
present, all solutions were evaluated, and a lower confidence flag 
was assigned (in almost all cases, such sources were classified as 
potential Ha emitters). 

We obtain redshifts for all NBj emitters with a usable spec¬ 
trum above a flux of 1.5 x 10“'^®ergs“^ cm“^, but we also ob¬ 
tain redshifts for more than half of sources with fluxes down to 
1.0 X 10“^® ergs“^ cm“®. In total, out of the 241 usable spec¬ 
tra, we detect a strong emission line and obtain redshifts for 185 
sources. For 135 sources, we are able to obtain a robust spectro¬ 
scopic redshift, based on at least 2 or more emission lines. Sources 
without a redshift determination are all lower flux sources, requir¬ 


ing higher exposure times in order to be significantly detected (note 
that we integrated for only 1.8ks per source). We show examples of 
FMOS spectra in Figure]^ which shows examples of Ha and [Out] 
emitters. 


3.4.2 MOSFIRE observations 

MOSFIRE observations were obtained on 5 November 2014, un¬ 
der clear conditions. The seeing was 0.8”. We observed two masks, 
both with 0.7” slits. The first one was centred on 22:13:24.89, -00 
12 16.31 (pi), while the second one was centred on 22 13 30.27, 
-00 53 21.24 (p3). Individual exposures were 120 s in each A and 
B positions, and we repeated each twice, for a total exposure time 
of 480 s per pixel for each of the masks. We reduced the data us- 
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Spectroscopic Redshift (z) 


Figure 5. Compaiison between photometric and spectroscopic redshifts for 
NB j selected line emitters for our S A22 NB survey. Redshifts for the emis¬ 
sion lines that we select are shown as vertical lines. Many other lines are 
not labelled, but found, e.g. from low to high redshift: Pa7l0870, NI10406, 
CI9853/9827, [SIII]9533, HeII8237, 017774 0117330/7320, [ArIII]7135, 
[SII]6731/6716, H«6563, [OI]6363/6300, HeI5876, [SIII]631 l-HHeII6311, 
[Olll]5007/4959, H/34861,[Oll]3727. 


ing the MOSFIRE team data reduction pipeline (DRlf^. The DRP 
produces flat-fielded and wavelength calibrated combined 2D spec¬ 
tra for the individual objects. The spectra are wavelength calibrated 
using the sky lines and then. In practice DRP follows very similar 
steps to those used to reduce FMOS data. 

For mask pi we were able to observe 17 of our candidate line 
emitters, obtaining robust redshifts for all 17 sources. For mask p3 
we targeted 8 of our line emitter candidates, and obtained a robust 
spectroscopic redshift for 7 out of the 8. For the missing source the 
S/N was too low to detect the emission lines. We note that for most 
sources the S/N obtained with MOSFIRE are at least comparable, 
and in most cases much higher than FMOS, despite MOSFIRE data 
having an exposure time of only ~ 13 per cent that of EMOS. How¬ 
ever, FMOS has a larger multiplexing, allowing to target about 7- 
10 times more sources (taking into account cross beam switching), 
and thus it is competitive with MOSFIRE for our targets, i.e., for 
sources spread over relatively wide areas and that have emission 
lines which avoid strong OH lines. 

In total, we obtained redshifts for 24 sources with MOSFIRE. 
We show examples of MOSFIRE spectra in Figure Figure 
shows the full spectroscopic redshift distribution. 

Together, the FMOS and MOSFIRE data-sets significantly 
add to the number of spectroscopic redshifts in the sample of line 
emitters, not only allowing for a spectroscopically confirmed sam¬ 
ple of just over 300 Ha emitters at 2 = 0.81, but also allowing 
to directly investigate the contamination by the adjacent [Nil] line 
and contributing to a robust correction (see j ]4.1^ . 

Finally, by compiling all the 511 spectroscopic redshifts, we 
test the accuracy of our photometric redshifts for line emitters, 
which can be seen in Figure]^ Figureshows that the redshifts 
overall agree well, but that there is some scatter, particularly for 
z > 1. We use our results, and our full spectroscopic sample, to 
optimise our selection of different line emitters and to estimate the 
completeness and contamination of each sample. 


https://code.google.eom/p/mosfire/ 


3.5 Testing SA22 photometric redshift selection using 
COSMOS 


In order to further check the validity/accuracy of using our EAZY 
photometric redshifts (pz) based on ugrizJK and our NB for the 
bulk of our sample, we make use of similar data available in the 


COSMOS field jSobral et al.|20T3a| i. As the COSMOS field ( |Scov-| 
|ille et al.|2007||Capak et al.|2007| McCracken et al.|2012| has been 
widely studied, there is data available in 30 bands (see e.g. |IIb^ 
|et al.|200^|2010[[2013^ , making photometric redshifts much more 
reliable, while there are also relatively more spectroscopic redshifts 
available in the literature. Thus, in order to have a fully compa¬ 
rable sample, we compute photometric redshifts in COSMOS by 
following the same method as for our SA22 sample, i.e., by using 
ugrizJK and NBj. This allows us to directly compare photomet¬ 
ric redshifts for a sample of similar line emitters and investigate any 
biases/incompleteness due to the limited availability of photometric 
bands. 


For the 700 sources selected as line-emitters in COSMOS by 
ISobral et al.|(|201 3a^ in their NB j band, 76 had spectroscopic red- 
shifts ( [Lilly et al.|2009) . As a first step, we derive photometric red¬ 
shifts using ugrizJK in COSMOS without the inclusion of the 
NBj filter and compare them with those derived with the inclu¬ 
sion of the NB. We find that both photometric redshift sets are in 
good agreement for most redshifts, except for the redshifts of strong 
emission lines such as Ha at 2 = 0.8; for these photometric red¬ 
shifts clearly benefit from the addition of the NB and are always 
more accurate at recovering spectroscopic redshifts, resulting in 
typical photometric redshift uncertainties of « 0.05. We compare 
the photometric redshifts obtained in COSMOS, with all bands, 
from |Ilbert et ^ j2009^ , with those we derive with EAZY using 
the restricted set of bands including NB. Our results are shown in 
Figure We find good agreement between both. We find that for 
the bright {NBj < 21.2, AB) subset of the sample, the EAZY 
photometric redshifts are consistent to « 0.02. For fainter sources, 
the EAZY photometric redshifts can sometimes differ significantly 
from those using all bands, due to larger error-bars for the faint 
sources, but also because the lack of bands redder than K and the 
lack of H and UV bands. 


Finally, using the data and computed photometric redshifts in 
COSMOS we check our selection of emitters that we will apply 
for SA22. We apply the same selection criteria (see §3.6| > in COS¬ 
MOS, using both the |Ilbert et al.| ( [200^ photo- 2 s based on 30-bands 
and our photometric redshifts with ugrizJK+^B. We find that 
our SA22-like selection criteria (photo -2 with just ugrizJK+NB 
and colour-colour selections; see Section [J!^ applied to COSMOS 
recovers 90% of the Ha emitters (97 % of the spectroscopically 
confirmed Ha emitters) that are found with much more accurate 
photometric redshifts. Our photometric redshifts tend to sometimes 
favour higher redshift solutions compared to [Ilbert et"^ ( |2009| l, 
which leads to 1 spectroscopic Ha emitter being selected as [Oil] 
emitter (but this is a small effect). However, apart from this, the 
selection is very clean and leads to the conclusion that for our 
emitters, ugrizJK+NB photometric redshifts -l- colour-colour se¬ 
lections are sufficient and comparable to the best photometric red¬ 
shifts in COSMOS. We check that similar conclusions are reached 
for H/3-l-[OlIl] and [Oil] emitters. Therefore, our samples selected 
from SA22 and COSMOS (all selected with the same selection cri¬ 
teria) are unlikely to have strong differences and should have very 
comparable completeness (~ 90%) and contamination (~ 10 %) 
fractions. 
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Figure 6. A comparison between photometric redshifts for the COSMOS 
sample of NBj line emitters (from [^bral et aLj^lOlBa) ), derived by us 
with EAZY using ugrizJK and NB magnitudes and those derived by|n] 
|bert et ar]j2QQ9^ with 30 bands from FUV to IRAC 8 fim. We find a good 
correlation between the two estimates, with a relatively small scatter be¬ 
low z < 1. The scatter increases at higher redshift, particularly due to the 
lack of IRAC bands, but our photometric redshift ranges used for selection, 
together with our colour-colour selections are able to tackle that. 


Figure 1 . B — Rvsi — K colour-colour separation jSobral et al.|2Q09^ that 
we use (together with photometric redshifts) to identify Ha emitters within 
the full sample emitters. Here, only narrow-band line emitters are plotted. 
Note how the selection is able to distinguish between 2 ~ 0.8 sources, 
those at 2: > 1 and those at lower redshifts (2 < 0.6), but also that only 
~ 5 % of all spectroscopically confirmed Ha emitters do not fall on our 
BRiK colour-colour selection, thus implying high completeness. 


3.6 Selection of Ha, [Oiii]/H^ & [Oil] emitters 

The selection of Ha, [Olll]/H/3 & [Oil] emitters at z = 
0.81,1.4, 2.2, respectively, is done following |Sobral et'^j2013^ , 
by using a combination of photometric redshifts (and spectro¬ 
scopic redshifts, when available) and colour-colour selections op¬ 
timised for star-forming galaxies at the redshifts of interest {z ~ 
0.8,1.4, 2.2). We apply these selections not only to the sample of 
emitters in SA22 presented here, but also to COSMOS and UDS 
(from [Sobral et al.|2013a[ t, in order to obtain larger samples over a 
larger number of independent volumes. Spectroscopic redshifts are 
also used to evaluate the completeness and contamination of the 
sample, although due to the selection function of most of the liter¬ 
ature spectroscopic redshifts they are only really useful for sources 
up to z ~ 0.8. The selection criteria applied here are the same as 
ISobral et al.H2013^ for Ha at z ~ 0.8, while for [OlIl]-l-H/3 at 
2 ~ 1.4 we use their criteria to select Hq sources at 2 ~ 1.47, 
and for [Oil] emitters at 2 ~ 2.2 we use their criteria for 2 ~ 2.2 
Hq emitters. We check with COSMOS and UDS that these criteria 
work well (see |Khostovan et al.|2015| for more details). 

The following Sections describe the specific selection criteria 
used to identify the different line emitters, hut briefly: i) we apply 
a photometric redshift range, centred on the expected redshift of 
the line, and that takes into account typical errors (motivated hy the 
estimated uncertainties in the photometric redshifts and on e.g. Fig¬ 
ures |^and|^ and the proximity of any other strong emission line, 
ii) we use colour-colour selections (see e.g. Figure |7]), to increase 
completeness and hi) we explicitly remove any spectroscopically 
confirmed contaminant and include any spectroscopic confirmed 
source in the sample. Finally, if a source is classified as e.g. Hq it 
can no longer he selected as a higher redshift emitter and, if clas¬ 
sified as [OlIl]-l-H/? it will not be able to be selected as an [Oil] 
emitter. We use our spectroscopic redshifts in order to test different 
photometric redshift cuts and study the completeness and contam¬ 
ination of our samples. We find that we need to use a wider cut in 
photometric redshifts than [Sobral et al.| j2013a|) in order to max¬ 


imise the completeness (95% for Ha emitters), while maintaining 
the contamination at a low level of ~ 10 — 15%. This a simple 
consequence of our photometric redshifts in SA22 having a larger 
uncertainty. 


3.6.1 Ha emitters at z = 0.81 


Hq emitters at 2 ~ 0.8 can he distinguished from lower and higher 
redshift emitters very efficiently by using i) our photometric red¬ 
shifts and ii) the BRiK ( [Sohral et al.|2009^ colour-colour diagram 


(see Figure in comhination with photometric redshifts and 


spectroscopic redshifts ( [Sohral et al.|2009al >. We classify a source 
as Hq if i) its photometric redshift is within Zphot = 0.70 — 0.95 
(1523 sources) or ii) if it satisfies the BRiK colour-colour criteria 
(Figure |7| and does not have 1.3 < Zphot < 1.5 (likely [Olll]) 
or 2.0 < Zphot < 2.4 (likely [Oil]): this step adds 1330 extra 
sources (a further 200 sources that satisfy the BRiK selection are 
not included due to their photo-zs clearly pointing towards higher 
redshift line emitters). We use the BRiK colour-colour selection 
to guarantee high completeness even for the faintest sources, where 
photometric redshifts become unreliable and/or unavailable. We re¬ 
cover 286 spectroscopically confirmed Hq emitters (see Figure 
for a few typical examples), remove 34 sources spectroscopically 
confirmed to be other emission lines either than Hq ([Sii] emit¬ 
ters are clearly the largest contaminant, but there are also a few 
low redshift emission lines for which VVDS and VIPERS is ex¬ 
tremely complete), and introduce 15 spectroscopically confirmed 
Hq emitters which did not satisfy the photometric selection cri¬ 
teria. Comparing this number with the 286 emitters spectroscopi¬ 
cally confirmed to be Hq which our selection successfully recov¬ 
ers indicates ~95% completeness. With the current spectroscopic 


For SA22 we estimate B band photometry by using g and mimicking 
the necessary empirical coiTection extracted from the COSMOS and UDS 
[Sobral et al.|2013^ , for which u,g and B data is available. The correction 
is: (B - r)oor = {g - r) X 1.1657 - 0.223 
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follow-up, we estimate contamination on the level of ~ 10%, being 
dominated by [Sll] emitters at a = 0.76. We note that this is a con¬ 
sequence of photometric redshifts and colour-colour selections not 
being able to completely distinguish between Ha and [Sll] emitters 
as they are very close in redshift. 

In total, we select 2834 Ha emitters in SA22 (photometric 
redshift distribution of the final sample is shown in Figure]^. Out of 
all Ha emitters in SA22, 295 are spectroscopically confirmed (see 
Figure|^for examples from FMOS, and Figurej^for the composite 
spectra), from VVDS (75), VIPERS (74), KMOS (20; [Sobral et'af] 
|2013b| [Stott et al.||2014|l, 10 sources followed-up with NTT and 
WHT, 16 with MOSFIRE and 100 with FMOS. 


3.6.2 [OIII] and HP emitters at z ~ 1.4 

In order to select [OlIl]-l-H/3 emitters^ we use the same selection 
as for Ha emitters at z = 1.47 in [Sobral et al.[p013a^ . Our se¬ 
lection criteria are the following: i) we disregarding all candidate 
Ha emitters obtained in the previous subsection (rejection of 2834 
sources), ii) we use our photometric redshifts and apply the selec¬ 
tion 1.0 < Zphot < 1.8 (911 sources) and iii) we use our BzK se¬ 
lection to identify z ~ 1.4 sources from which we remove sources 
likely to lie at 2 > 2 on the basis of their izK colours (see [So-[ 
[bral et al.|2013a^ and a further 41 sources with 2.0 < Zphot <2.4 
which are likely to be [Oil] emitters; this adds 141 sources to the 
sample. Just like for Ha, we also use spectroscopic redshift infor¬ 
mation to remove eight non-[OlIl]-l-H/3 emitters in the sample and 
add twelve spectroscopically confirmed [Oin]-l-H/? emitters. Based 
on the limited spectroscopic information, we estimate ~ 80% com¬ 
pleteness and ~ 15% contamination. We have a total of 46 spectro¬ 
scopically confirmed sources; in Figure|^we show some examples. 
By following these selection criteria, we select 1056 [OlIl]-l-H/3 
emitters in SA22 (photometric redshift distribution of the final sam¬ 
ple is shown in Figure]^. We note that our limited spectroscopic 
follow-up already allows us to constrain the fraction of H/3 and 
[Olll] emitters within this sample. Within the full sample of spec¬ 
troscopically confirmed [OlIl]-tH/3 emitters, we find that ~ 16 per 
cent are H/1 emitters (HP line detected by the NB filter), and that 
H/3 emitters have lower luminosities than [Olll] emitters, likely 
contributing more towards the faint-end of the luminosity function, 
but very little at the bright end. We discuss the implications in more 
detail in §4.2[ 


3.6.3 [Oil] emitters at z = 2.2 

We select [Oil] emitters at 2 = 2.2 with a similar selection as 
for Ha emitters at 2 = 2.23 in [Sobral et al.[ ( |2013a[ >. Our selec¬ 
tion is as follows: i) we disregard all candidate Ha emitters and 
[OlIl]-l-H/3 emitters obtained in the previous subsections (rejection 
of 3888 sources); ii) we use our photometric redshifts, applying the 
selection 1.8 < Zphot < 2.7 (424 sources); and iii) we use a BzK 
cut, to isolate 2 > 1.5 galaxies jSobral et al.[2013al >, adding 69 new 
sources to the sample. We also use spectroscopic redshift informa¬ 
tion to remove seven non-[Oll] emitters in the sample and introduce 
two spectroscopically confirmed sources. We obtain a final sample 


It is not possible to completely distinguish, with photometric redshifts 
and/or colour-colour selections between [OlII] and H/3 emitters due to a 
very similar redshift, although we expect [OlII] emitters to dominate the 
sample - see 34.2[ 


Table 3. Summary of the information in the SA22 catalog. It first shows the 
number of NB j detections, then the number of sources selected as emitters 
in Section 3.1. It shows the number of candidate emitters after removing 
spurious sources, the number of candidate stars and the final number of ro¬ 
bust emitters. The remaining columns present the number of Hq, [OlII]/H/3 
and [On] candidates. When available, the number of spectroscopic redshifts 
are shown as well. 


Sample 

No. of sources 

2 -spec 

NBj detections 

346,244 

16964 

Candidate emitters (after visual check) 

6315 

541 

Stars 

339 

30 

Robust emitters 

5976 

511 

Ho (2 = 0.81) 

2834 

295 

[Oni]+H/3 (2 = 1.37/1.44) 

1056 

46 

[On] (2 = 2.18) 

488 

19 

[Sll] (2 = 0.76) 

- 

37 

-2 < 0.7 ( 2 : > 3 or unidentified) 

1638 

82 (26) 


of 488 [Oil] emitters in SA22, with 19 sources being spectroscopi¬ 
cally confirmed. Our spectroscopic sample is extremely limited for 
our [Oil] sample, but based on the available redshifts, we estimate a 
~ 90% completeness and a < 25% contamination. We note that we 
do not apply the Lyman-break selection in order to exclude 2 > 3 
emitters as in [Sobral et al.[p013a[ l, because no major emission lines 
are found at slightly higher redshift — but the reader is referred to 
[Matthee et al.[ ( [20f4l > for a full discussion of emitters at much higher 
redshift. 


3.6.4 Final Samples: SA22 

Using the data and steps mentioned in previous chapters, we 
add an identifier in the range 0 to 3 to classify between other 
lines/unclassified (0), Ha (1), [OlIl]-l-H/3 (2) and [Oil] (3) to 
the catalogue of emitters containing the final 5976 robust excess 
sources. Table [^presents the final numbers per sample, including 
the number of spectroscopically confirmed. 


3.6.5 Samples in COSMOS and UDS 


We supplement our sample using the HiZELS catalogue of line 
emitters available from [Sobral et al.| ( [2013a^ , which has been de¬ 
rived in a very similar way to ours, and with almost identical data 
quality and d epth. For Ha, we use the samples presented by [Sobra]] 


et al. 


1 2Q13a| >: 425 Ha emitters in COSMOS and 212 in UD^f 


In order to select [OlIl]-bH/3 and [Oil] emitters, we apply 
exactly the same criteria as in SA22. In COSMOS, we find 159 
[OlIl]-bH,fl and 41 [Oil] emitters. In UDS, we find 128 [OlIl]-l-H,5 
emitters at 2 = 1.4 and 43 [Oil] emitters at 2 = 2.2. We note that 
the selection criteria (optimised for the COSMOS and UDS fields) 
in [Khostovan et al.| ( [2015t to select [OlIl]-l-H/3 and [Oil] emitters 
with the HIZELS NB,/ sample is slightly different than ours in the 


By applying exactly the same selection as in SA22 we would obtain two 
extra sources in COSMOS, and three extra sources (and less six) in UDS 
relative t o[Sobral et al.[ j2013a] , but we chose to use the sample presented in 
[Sobral et al.[j2013a] for consistency. These minor differences in the sample 
make no difference in the results, particularly as they are found at the fainter 
end. 
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Table 4. A summary of our final samples of Ho, [0lII]+H/3 and [Oil] emit¬ 
ters at z ~ 0.8, z ~ 1.4 and z ~ 2.2, respectively. 


Samples 

No. of sources 

Volume 

Depth (logioL) 

Ha {z ~ 0.8) 

(#) 

(10® Mpc®) 

(ergs-l) 

SA22 

2834 

0.9 

41.35 

COSMOS 

425 

0.1 

41.25 

UDS 

212 

0.08 

41.30 

Full Sample 

3471 

1.1 

41.35 

[Olll]-(-H/5 (2 ~ 1.4) 

(#) 

(10® Mpc®) 

(ergs“l) 

SA22 

1056 

2.48 

42.0 

COSMOS 

159 

0.33 

41.90 

UDS 

128 

0.24 

41.95 

Full Sample 

1343 

3.05 

42.0 

[On] (2 ~ 2.2) 

(#) 

(10® Mpc®) 

(ergs“l) 

SA22 

488 

2.62 

42.6 

COSMOS 

41 

0.29 

42.50 

UDS 

43 

0.22 

42.55 

Full Sample 

572 

3.12 

42.6 


photo-z selection (slightly more restrictive due to photo-zs in COS¬ 
MOS and UDS being better than in SA22), but although it results 
in some very minor differences in the source numbers (maximum 
~ 15% but typically within 5-10%), the samples are essentially the 
same (> 90% of the samples are the same) and lead to the same re¬ 
sults. Here we chose to use the same selection criteria as for SA22 
to have fully consistent samples across fields, but we note that even 
if we applied [Khostovan et al.|p015| l selection criteria our results 
would not change. 

We also note that the HiZELS data jSobral et al.|2009^|2013a^ 
are slightly deeper (by ~ 0.1 dex in luminosity) than our SA22 
survey and that the narrow-band filter used by |Sobral et al.|p013a^ 
is slightly wider, and thus naturally recovers a higher number of 
emitters per deg^ if that is not taken into account. 

We provide a summary of the sample in the three fields in Ta¬ 
ble 13 The final samples (SA22, COSMOS and UDS) are by far 
the largest ever assembled, yielding 3471 Ha emitters aX z ~ 0.8 
(~ 400 spectroscopically confirmed), 1343 [OlIl]-l-H/3 at a ~ 1.4 
and 572 [Oil] emitters at z ~ 2.2. Some relatively significant vari¬ 
ations in source densities are found across fields and within fields. 
This is evaluated in §5, where we present how important is cosmic 
variance on different scales for each one of the emission-lines. 


luminosities (so they are not the highest luminosity line emitters, 
but none is a weak line emitter). 

We also use deep IRAC data, in both COSMOS and UDS, to 
look for sources dominated by signatures of the stellar bump (star¬ 
forming dominated), thus showing a blue colour beyond ~ 1.6 /rm 
rest-frame and sources dominated by signatures of a red power-law 
(AGN dominated), showing red colours even beyond ~ 1.6/rm 
rest-frame (see e.g. |Garn et al.|20T0| ). Given the typical errors on 
IRAC photometry, particularly for redder bands, and motivated by 
the colours of our X-ray AGN (in order to recover the majority 
of them with the appropriate IRAC colours given the redshift of 
the sources), we apply the following cuts to select potential AGN, 
which should be mostly interpreted as an upper limit to the AGN 
contamination. For Ha {z = 0.8): [3.6]-[4.5]> 0.0 (10% potential 
AGN); for [OlIl]-l-H/3 (z ~ 1.4): [4.5]-[5.8]> 0.1 (18% poten¬ 
tial AGN) and for [Oil] (z = 2.2): [5.8]-[8.0]> 0.2 (23% poten¬ 
tial AGN). The contamination for these relatively “typical” sources 
can be well approximated by a constant contamination across lu¬ 
minosities at least up for the luminosities probed by the UDS and 
COSMOS samples. 

In order to constrain the fraction of AGN among the most lu¬ 
minous Ha emitters, we use the results from |Sobral et al.| ( |2015| l. 
ISobral et al.| ( |2015| l present the results of the spectroscopic follow¬ 
up of very luminous Ha emitters, allowing to constrain the fraction 
of AGN for Ha luminosities L > L* (up to L ~ 50L*). Their 
results show that while the AGN fraction is relatively constant up 
to L ~ L* (around 10%), it strongly correlates with L/L* for 
higher luminosities, reaching ~ 100% by L ~ SOL*. Therefore, 
for Ha emitters, we use results from X-rays and from IRAC, which 
allow us to have a global AGN fraction for typical luminosities, and 
the results from |Sobral et al.| ( |2015t , which allow us to estimate the 
AGN fraction at the highest luminosities (for Ha emitters). We thus 
assume that our Ha sample {z = 0.8) will likely be ~ 10% con¬ 
taminated by AGN up to L ~ L*, and we use the AGN fraction as 
a function of L/L* found by|Sobral et al.j(|2015), which is given 
by 0.59 X logio(L/L*) + 0.112. 

For our sample of [OlIl]-tH/3 and [Oil] emitters we have no 
information regarding the AGN fraction of the highest luminosity 
emitters. We thus assume a constant AGN fraction when convert¬ 
ing luminosity densities to star formation rate densities. For our 
sample of [OlIl]-l-H/3 we assume we will be ~ 15 % contaminated 
by AGNs ( [Stott et al.|2013b[ find that for Ha emitters at 2 : = 1.4 
only 10% are AGN) and for our [Oil] sample at 2 = 2.2 we use a 
correction of ~ 20 %. 


3.7 AGN fraction for line emitters 


We take advantage of the similar samples we selected in COS¬ 
MOS and UDS to estimate the fraction of potential AGNs in our 
SA22 sample (where no Chandra data and no IRAC data are avail¬ 
able). We start by using C-COSMOS ([Elvis et al.||2009) to find 
that, in agreement with [Garn et ^ j2010^ , only ^ 1% of Ha 
sources are detected in the X-rays (5 sources. X-ray luminosities of 
^q 42 . 8 ±o.ii gj.g ^- 1 ^ Pjjj. [Oiii]-hH; 9 and [Oil] samples in COS¬ 
MOS, we only find one X-ray match per sample (both sources with 
X-ray luminosities of ~ erg s“^), thus resulting in a very 

small fraction (< 1 % for [OlIl]-(-H/3 at 2 ~ 1.4) and ~ 2 — 3% 
for [On] at 2 = 2.23. It is therefore clear that the fraction of X-ray 
detected AGNs in our samples are only at the level of ~ 1 %. All 
X-ray AGN within our sample are found to have intermediate line 


4 METHODS: COMPLETENESS AND CORRECTIONS 

4.1 Flux corrections: [NII]/Ha correction 

When computing line fluxes and equivalent widths for the Ha emit¬ 
ters (Ha line), one must note that the adjacent [Nil] lines at 6548A 
and 65 83A will also contribute to both quantities, increasing them 
both (see e.g. Figurej^l. [Nll]6583A is clearly the strongest and the 
most important to take into account (see Figurej^l; we refer to it as 
simply [Nil] for the remaining of the paper. One way to correct for 
this is to use the SDSS relation between F[nii]/Fhci and the total 
measured equivalent width EW(Ha-(-[Nll]), ( [Villar et al.[2008[[So-[ 
[bral et al.|2012^ . Given the representative spectroscopic follow-up 
of our Ha emitters, we can directly test whether the polynomial 
correction presented in [Sobral et al.[p0T^ (based on SDSS) is ap¬ 
propriate for 2 ~ 1. 
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Restframe Wavelength (A) 


Figure 8. Stacked spectra of 2 = 0.81 Hq emitters obtained with 
FMOS and MOSFIRE. As a whole, our Ha emitters have a metallic- 
ity 12 + logio(0/H) = 8.56 di 0.05, slightly sub-solar (but consis¬ 
tent with solar metallicity 12 + logio(0/H) = 8.66 di 0.05). We 
also make significant detections of the [SlI] doublet, with a line ratio of 
I([Sii] 6716 )/I([Sii] 673 i) = 1-33 di 0.08, implying an electron density of 
40-200cm“^. The [SlI] 67 i 6 /I(Ha)= 0.14 di 0.02 ratio also implies ioni¬ 
sation potential of logio(U) = —3.9 di 0.5 cm“^. 



Rest-frame EWho+inii] (A) 

Figure 9. The observed anti-correlation between the [NlI]6583A/Ho line 
ratio and the rest-frame EW(Ha-l-[NlI]6583A) at 2 : = 0.81 in this study 
and a comparison with the polynomial fit to SDSS at low redshift. Given 
the range of EWs and line ratios probed, we find that while the SDSS fit is a 
good fit to the data at z = 0.81, a linear fit does equally well, and is a much 
simpler way to robustly estimate [NII]/Ha. The relation is /([NII]/Ha) = 
—0.296 X logxo(EWHa-|-[Nll]) + 0.8. We show three bins which split the 
sample in 3 relatively evenly in parameter space and that do not overlap, 
together with the stack for the full sample, and 3 other stacks which split 
the sample in terms of number of Ha emitters. 


We use our FMOS and MOSFIRE data and stack as a function 
of rest-frame EW(Ha-l-[Nll]) (measured from the NB data as that 
is what we want to evaluate). We recover a clear anti-correlation 
between F[Nn]/FHa and rest-frame EW(Ha-l-[Nll]), as shown in 
Figure]^ We find that the trend is fully consistent with SDSS of 
decreasing F[Nii]/FHa as a function of EW(Hq:-i-[Nii]). We note, 
as shown in Figure]^ that a linear relation is an even simpler cor¬ 
rection which is extremely similar to the SDSS relation. The linear 


relation that we derive is: 


/([NII]/Ha) — —0.296 x logj^g(EWHa-i-[Nii]) -1-0.8 (8) 

and is valid for rest-frame EW(Ha-l-[Nll]) from ~ 15 A to ~ 600 A 
(see Figure]^. 

We note that if a single correction is applied (for similar sur¬ 
veys), then the value to be applied for FfNiij/FHc should be 25%, 
as the FfNii]/FHa for the full stack (median stack, see Figure[8| is 
0.25 ± 0.05. This suggests (in agreement with e.g.|Swinbank et aE 


2012a[[Stott et al.|2013b[ [Sobral et al.|2013b| l, using the Pettini & 


Pagel ([2004 1 calibration, that as a whole, our Hq emitters have a 


metallicity 12 -|- logic(0/H) = 8.56 ± 0.05, slightly sub-solar 
(solar metallicity 12 -|- logic (0/H) = 8.66 ± 0.05). 

In our stacked spectra, we also make significant detec¬ 
tions of the [Sll] doublet. We find a median line ratio of 



we emphasise that the stacked spectra provide only information on 
the median value, and no indication on the range of values. 


4.2 The relative contributions from [OiiiIscct, [Oiii]4959 and 
FI/3 to the sample of [Oiii]-tH/3 emitters at 2 ~ 1.4 

For [OlIl]-bH/3, we start by noting that based on our spectroscopic 
follow up, ~ 16% of the [OlIl]-(-H/3 emitters turned out to be H/3 
and the rest to be [Olll]. We also find that spectroscopically con¬ 
firmed H/J emitters are found to have lower luminosities than spec¬ 
troscopically confirmed [OlIl] emitters. This, combined with the 
low fraction within the sample, would be a relatively good motiva¬ 
tion towards ignoring such emitters (this is done by most studies, 
e.g. [Ly et al.|[2007[ ). However, attempting to further split the H/3 
and [Out] luminosity functions, we take the following approach. 
We use the [Sobral et al.[j2013a[ l 2 = 1.47 Ha luminosity func¬ 
tion, and, with a simple assumption of Aho = 1 mag and case B 
recombination, we predict the H/3 Luminosity function at 2 ~ 1.44 
(also taking into account what flux we would recover with our NB 
filter). We find that the number densities of H/3 emitters are signif¬ 
icantly below the number densities of [OlIl]-(-H/3 emitters. Given 
our depth, this simple prediction is very much in line with the 16% 
fraction of H/3 found (and with H/3 emitters being preferentially 
found at fainter luminosities). 

Due to the separation in wavelength, it is not possible for both 
H/3 and [Olll] to contribute to the NB flux measured. We note, 
nonetheless, that there is a narrow range of redshifts where both 
lines from the [Out] doublet are detected (see one example of an 
emitter in Figurej^where both [Olll] lines contribute to the narrow- 
band flux) at opposite wings of the filter. Indeed, the [OlIl]4959 
emission line is found to be even more important within the sample 
than H/3 and with important consequences for deriving luminos¬ 
ity functions. While we find that [Olll]5007 (only) dominates the 
sample of spectroscopically confirmed [OlIl]-l-H/3 emitfers, repre¬ 
senting ~ 50% of the full sample, [OlIl]4959 (only) represents 
27% (see e.g. Figure]^, almost twice as common as H/3 and 7% 
of the sources are actually detections of both [Olll] lines at op¬ 
posite wings of the filter (e.g. Figure]^. Within the spectroscopic 
sample of [OlIl]-l-H/3 emitters, we study the fraction of each emit¬ 
ter as a function of luminosity and show the results in Figure [T^ 
Our analysis of the spectroscopically confirmed emitters shows that 
[Olll]5007 emitters are distributed in a very Schechter like form, 
at all luminosities, H/3 are only found at lower luminosities, but 
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Figure 10. The fraction (within the spectroscopic sample) of [OlII]5007, 
[OlII]4959 and H/3 emitters as a function of [OlII]+H/3 luminosity within 
the [OlII]+H/3 sample of emitters at 2 : ~ 1.4. We find that [OlII]5007 emit¬ 
ters dominate at all luminosities, and that H/5 emitters only start to have a 
small contribution for lower luminosities. 
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density. We also note that for the [Oil] sample, we do not need to 
apply or investigate any correction, as no other line (close enough 
in rest-frame wavelength) is expected to contribute to the flux mea¬ 
sured. 


4.4 Completeness corrections 

Fainter sources and those with weak emission lines might be missed 
and thus not included in the sample; this will result in the underesti¬ 
mation of the number of emitters, especially at lower luminosities. 
In order to account for that we follow |Sobral et ar] ( |2013a^ to esti¬ 
mate completeness corrections per sub-field per emission line. Very 
briefly, we use sources which have not been selected as line emit¬ 
ters (E < 3 or EW< 30 A) and that have a photo-a within ±0.5 
of the appropriate line of interest (thus, we use different samples to 
estimate the completeness of different lines; see e.g. |Sobral et al.| 
l |2012^ ). We then add emission-line flux to all those sources, and 
study the recovery fraction as a function of input flux. We do these 
simulations in a sub-field by sub-field basis. We then apply those 
corrections in order to obtain our completeness-corrected luminos¬ 
ity functions. 


[OlIl]4959 emitters, and simultaneous detections of [Olll]5007 and 
[OlIl]4959 lines, are found both at the faintest luminosities (in the 
NB) and at the brightest luminosities, thus significantly boosting 
number counts at both the faintest luminosities (likely steepening 
the faint-end slope) and at the brightest luminosities (likely con¬ 
tributing to a non-Schechter form). 

Our results are very important in order to interpret the results 
from the luminosity function in [5.1.2 and for other similar stud¬ 
ies (e.g. |Ly et al.]|2007[ |Khostovan et al.|[2015| l. Thus, while we 
present a [OlIl]-l-H/3 luminosity function, we caution that in order 
to properly derive luminosity functions for each of the individual 
lines contributing to it, a much more detailed spectroscopic follow¬ 
up is needed, and that simply assuming that similar samples will be 
dominated by [Olll]5007 (ignoring, for example, the [OlIl]4959 
line, or simultaneous measurements of the doublet), will lead to 
strong systematic errors/biases. We therefore refer to these emitters 
as [OlIl]-l-H/3 throughout the paper. We also take our results into 
account when estimating the total volume probed by our survey for 
[OlIl]-l-H/? emitters. Given that [Olll]5007 dominates the sample, 
we use the full volume probed with that emission line, but we cor¬ 
rect for the expected addition of extra emitters (picked up over extra 
volumes). Thus, we also add to the total volume probed (because 
we are also sensitive to the other lines, and following our spectro¬ 
scopic results) 16% of full H/3 volume probed by us and 25% of 
the volume probed for [OlIl]4959 emitters. 


4.5 Filter profile corrections 

The narrow-band filter transmission function is not a perfect top- 
hat, so the real volume surveyed is a function of intrinsic lumi¬ 
nosity. For example, luminous line emitters will be detectable over 
a larger volume than the fainter ones, as they can be detected in 
the wings of the filters (although they will be detected as fainter 
sources in these cases). Low luminosity sources, however, will only 
be detectable in the central regions of the filter, leading to a smaller 
effective volume. 

In order to correct for this when deriving the luminosity func¬ 
tions, we follow [Sobral et aL] ( |2012^ . Firstly, we compute the lumi¬ 
nosity function assuming a top-hat narrow-band filter. We then gen¬ 
erate a set of 10^° line emitters with a flux distribution given by the 
measured luminosity function, but spread evenly over the redshift 
range being studied (assuming no cosmic structure variation or evo¬ 
lution of the luminosity function over this narrow redshift range). 
We fold the fake line emitters through the top-hat filter model to 
confirm that we recover the input luminosity function perfectly. 
Next, we fold the fake line emitters through the real narrow-band 
profiles - their measured flux is not only a function of their real 
flux, but also of the transmission of the narrow-band filter for their 
redshift. The simulations show that the number of brighter sources 
is underestimated relative to the fainter sources. A mean correction 
factor between the input luminosity function and the one recovered 
(as a function of luminosity) was then used to correct each bin. 


4.3 Luminosity Calculations and Extinction corrections 

In order to calculate luminosity functions for our samples, line 
fluxes are converted to luminosities by applying: 

Lline = 47rDLFline (9) 

where Dl is the luminosity distance. We use DL=5367Mpc for 
Ha emitters at z = 0.81, Dl= 9752.7 Mpc for [OlIl]/H/3 emitters 
at 2 = 1.4 and Dl= 17746.5 Mpc for [Oil] emitters at 2 = 2.2. 

We note that all our luminosity functions are observed lumi¬ 
nosity functions, not dust corrected. We only apply extinction cor¬ 
rections when converting luminosity density to star formation rate 


5 RESULTS 

5.1 Luminosity functions 

The estimate of the source density in a luminosity bin of width 
A (log L) centred on log Lc is given by the sum of the inverse vol¬ 
umes of all the sources in that bin, after correcting for complete¬ 
ness. The volume probed is calculated taking into account the sur¬ 
vey area and the narrow-band filter width, followed by applying the 
appropriate real filter profile corrections obtained in §4.5| 

The luminosity functions presented here are fitted with 
Schechter functions defined by the three parameters: a, (f)* and L* 




































14 D. Sobral et al. 


(see Table 1^. For the Ha luminosity function, we can still get a 
reasonable constraint on a, but the data are too shallow to do the 
same for [OlIl]+H/3 and [Oil], and thus for those we concentrate 
on just fitting for cfi* and L* (and fixing a to the values commonly 
used in the literature for comparison). The final luminosity func¬ 
tions are presented in Figures[TT]andp^ and in Tables [bT||B 2| and 
|B3[ We note that we do not apply any dust correction prior to fitting 
the luminosity functions - dust corrections are only applied when 
converting luminosity densities to estimates of the star formation 
rate density. The best fit Schechter functions are presented in Table 

a 


5.1.1 Ha luminosity function at z = 0.81 

We present our Ha luminosity function at 2 = 0.81 in Figure [TT] 
and Table [BT] We confirm a strong evolution in the Ha luminos¬ 
ity function from 2 = 0 to 2 = 0.8. We also find good agree¬ 
ment with both |Sobral et al.H2013a| l and |Ly et al.| ( |20lT| l, as can 
be seen by directly comparing the data points and the best fit lumi¬ 
nosity function from e.g. [Sobral et al.|f2013a[ l. However, as Figure 
|TT| clearly shows, we are able to not only extend the LF to higher 
luminosities, and to constrain the number densities of the most lu¬ 
minous Ha emitters with unprecedented accuracy, but also deter¬ 
mine it with 3x more bins while reducing the Poissonian errors 
significantly at every bin. While the agreement is relatively good 
at all luminosities, we find a slightly lower L* or slightly lower 
(f>*, in line with e.g. [Colbert et al.| ( (2013) . Our results in §5.3 also 
show that we over-come cosmic variance, with the errors becoming 
a few times smaller than all previous surveys, and that the differ¬ 
ences with respect to |Sobral et al.[p013a[ ) and [Ly et al.[^20lf| are 
fully explained by cosmic variance. Our results also confirm that 
the faint-end slope is relatively steep (we find a = —1.6 ± 0.2), 
in excellent agreement with both [Sobral et al.[ ( [2013a^ and [Ly et al.[ 
( |2011) and that L* evolves significantly from 2 ~ 0 to 2 ~ 0.8. 
We also note that our spectroscopically confirmed sub-sample (301 
spectroscopically confirmed Ha emitters in the SA22 sample) fully 
supports our results. Moreover, we find that the density of the most 
luminous emitters is still consistent with [Sobral et al.[ ( [2013a[ l, al¬ 
though, as shown in [Sobral et al.[ j2015^ , the most luminous Ha 
emitters at 2 ~ 0.8 have a significant fraction of AGN, particularly 
broad-line AGN. 


5.1.2 [OlllJ+HP luminosity function at z ~ 1.4 

We present the observed (non-dust corrected) [OlIl]-l-H/3 luminos¬ 
ity function at 2 ~ 1.4 in Figure [T^ and in Table [B^ We have 
fixed a = —1.6 for the best-fit. We find a clear luminosity evo¬ 
lution from 2~0to2~1.4 when compared to other luminosity 
functions from the literature, in line with what is seen for Ha. Inter¬ 
estingly, if one uses empirical conversions from [Olll] to Ha (e.g. 
[Ly et al.|2007[ l, [OlIl]/Ha « 1.05, and plot the Ha LF at 2 ~ 1.4, 
we find that it resembles the [OlIl]-l-H/3 luminosity function sur¬ 
prisingly well. In other words, [OlIl]-l-H/3 emitters at 2 ~ 1.4 have 
similar number densities to Ha emitters at the same redshift. How¬ 
ever, as noted in §4.2[ the [OlIl]-l-H/3 luminosity function is far from 
being simple/easy to interpret (as it is a composite of a few different 
lines, unlike Ha). 

Comparing with the largest HST slitless grism survey (|Col-| 
[bert et al.|2013^ , we find very good agreement at moderate to high 
luminosities, but we recover more line emitters at the faintest lu¬ 
minosities, as discussed earlier in §4.2[ Based on our spectroscopic 


data, our [OlIl]-l-H/3 luminosity function is likely a sum of three dif¬ 
ferent luminosity functions, all contributing in varying forms across 
the observed luminosities, thus showing a non-Schechter like form: 
[OlIl]4959, [OIII]5007, [0lll]4959-l-[0lll]5007 (where both lines 
are both being detected and measured by the narrow-band filter) 
and H/3. 


5.1.3 [Oil] luminosity function at z = 2.2 


We present the observed (non-dust corrected) [Oil] luminosity 
function at 2 = 2.2 in Figureand in Table [B^ This is the first 
time the [Oil] luminosity function has been measured at 2 = 2.2, 
and also constitutes by far the largest sample of [Oil] emitters at 
2 ~ 2. Since we do not constrain the faint-end slope, we fix it to 
the value measured at 2 ~ 1.5 by|Ly et al.[([2007[l; [Sobral et al.[ 
l [2m^ , a = -0.9. 


We compare our results with results from the literature at 
lower redshifts: [Galle g o et al.[j2002[ l, Takahashi et al. ^2007^ , [I^ 
[et al.[ ( [2007^ , |Sobral et al. ( 2012{ and Bayliss et al. j2011 1 . Similarly 
to what has been found with Ha and [OlIl]-l-H/3, we find a signif¬ 
icant L* evolution in the [Oil] luminosity function. However, for 
[On] we find an even stronger/faster L* evolution as L* becomes 
~ 100 X brighter from the local Universe to 2 = 2.2. We also 
show data from a recent study by [Khostovan et al.[ ( [2015[ l, based 
on the HiZELS (COSMOS-l-UDS) data. We find excellent agree¬ 
ment when comparing to [Khostovan et al^ ( [2015[ l, which probes to 
slightly lower luminosities, while we probe a much larger volume. 
Indeed, our fit to SA22 fits the results from [Khostovan et al.[j2015l > 
well. However, even when combined, our data-sets are still unable 
to constrain the faint-end slope, and thus it is possible that the faint- 
end slope is steeper than a = —0.9. 


5.2 Predictions for Euclid and WFIRST 


We find that approximately ~50 per cent of the emitters in our 
sample (for our full sample down to our flux limit) are likely Ha 
emitters at 2 = 0.81, while [OlIl]-fH/3 account for ~ 18 per cent 
and [Oil] for about 9 per cent. The remaining ~ 23 per cent con¬ 
sist of rarer emission lines (mostly at 2 < 0.8, but also at 2 > 3; 
see [Matthee et al.[p014[ l) and extra stars (with strong spectral fea- 
tures ). If our flux limit wa s similar to the planned Euclid wide sur¬ 
vey iLaureijs et al. |2012 |> of > 3 X 10 ergs we would ob¬ 
tain a lower fraction (but still significant) of line emitters which 
are not Ha, [OlIl]-l-H/3 or [Oil]: 18 %. Our Euclid-hke sample is 
dominated by Ha emitters (67%), followed by [OlIl]-l-H/3 emitters 
(12%), but still with some [Oil] emitters (3%). The raw (observed) 


(A ~ 1.2/rm) is 8.66x10“® Mpc“®, while at the same wave¬ 
length the observed number density of [OlIl]-l-H/3 emitters is about 
5 times lower (1.6x 10“® Mpc“®), and the number density of [Oil] 
emitters (0.3xl0“® Mpc“®) is almost 30 times lower than Ha at 
A ~ 1.2 fim. 

Using spectroscopic redshifts from the literature and from our 
own follow-up, we can have an even more robust quantification of 
the range of line emitters in the full sample. Figurej^shows the full 
distribution of spectroscopic redshifts. Figure[^also shows the dif¬ 
ferent spectroscopic redshift sub-samples obtained when restricting 
our sample to a much higher narrow-band colour significance (E). 
For E > 3 (our full sample), and based solely on spectroscopic 
redshifts, we find 59% Ha emitters, 9% [OlIl]-l-H/3 emitters, 7% 
[Sll] emitters, 4% [Oil] emitters, while 16% are other 2 < 1.0 
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Figure 11. Observed Ha luminosity function at z = 0.81 from our survey, comparison with previous surveys at a similar redshift, and evolution of the 
luminosity function (none of the luminosity functions have been corrected for dust extinction). We find a very good agreement between our results and those 
in the literature (well within the errors excepted from cosmic variance), but are able to extend the Ha luminosity function to much higher luminosities, as well 
as reducing the Poissonian errors significantly. The errors due to cosmic variance are < 10% at most. Our results agree with the strong L* evolution with 
increasing redshift. 



(including [Nil]) and the remaining 5% are other z > 1 emit¬ 
ters. Sub-samples with higher E will be samples with higher flux 
limits, thus changing the distribution of emitters within such sub¬ 
samples. For our sample, E > 4 corresponds to a cut in flux of 
1.3 X 10“^® ergs“^ cm“^, while the cut at E > 5 corresponds to 
cutting the flux down to a flux limit of 1.55 x 10“^® erg s“^ cm“^. 
For E > 4 (E > 5), based on spectroscopic redshift, we find 
74% (81%) Ha emitters, 9% (7%) [Oni]-tH/3 emitters, 6 % (4%) 
[Sll] emitters and 3% (3%) [Oil] emitters, while 6 % (3%) are other 
2 < 1.0 (including [Nil]) and the remaining 5% (2%) are other 
z > 1 emitters. It is therefore clear that a selection with a higher 
flux limit will be more and more dominated by Ha emitters, and 
thus our spectroscopic sample shows that both Euclid and WFIRST 
can expect ~ 80% of all line emitters to be Ha. 


5.3 Luminosity density and Cosmic star formation rate 
density 

By using the largest samples ever assembled and empirically com¬ 
puted uncertainties due to cosmic variance, we integrate our ob¬ 
served luminosity functions. We provide two different measure¬ 
ments for each emission line: i) the full analytical integral and ii) 
the numerical integration down to the observational limit. We pro¬ 
vide the observed luminosity density values in Table 

Prior to converting our observed luminosity densities for each 


line to a star formation rate density (psfr), we apply two cor¬ 
rections. First we do a simple correction for dust extinction. We 
assume 1 magnitude at Ha (Aho = 1). Dust corrections for 
[OlIl]-l-H/? and [Oil] are discussed below. Secondly, we assume a 
10-20% AGN contribution to our samples, following our results in 



In order to convert luminosity densities to star formation rate 
densities we use, for Ha (|Kennicutt|1998t: 


SFR(M 0 yr“^) = 7.9 x lO""*^ Lhc [ergs' 


( 10 ) 


For [OlIl]-l-H/3, as seen before, directly interpreting as a star- 
formation indicator has a significant number of caveats. However, 
cautioning the reader that any values need to be interpreted with 
caution and are mostly indicative, we use the relation derived from 
[Osterbrock & Ferh^ 1 2006|^jn good agreement with the empiri- 
cal calibration by|Ly et al.H2007|l, who found typical line ratios of 
~ 1 : 


SFR(M 0 yr = 7.35 x 10 L[oiii]+h/3 [ergs ^]. (11) 

For [Oil] we use the standard [Kennicutt|jl998^ calibration of [Oil] 
as a star formation tracer (calibrated using Ha): 

SFR(M 0 yr-') = 1.4 x lO”"' L[oii] [ergs”']. (12) 
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Figure 12. Left: Observed [OlII]+H/5 luminosity function at 2 ~ 1.4 from our SA22 survey and comparison with other studies/surveys (including the LF from 
UDS+COSMOS derived by |Khostovan et al.|2015) . Our results clearly show that the [OlII]+H/3 luminosity function is evolving very strongly with redshift, 
mostly due to a strong increase in L*, similarly to the behaviour of the Ho; luminosity function. Our best-fit is very similar to the observed Ho luminosity 
function at the same redshift (also shown), and, at the faint end, reveals a higher number of [OlII]+H/3 emitters than in |Colbert et al.|j2Q13t - these are likely 
explained by a combination of H/3 emitters (we also show our predicted H/3 luminosity function based on the Ha luminosity function from|Sobral et al.| 
j2Q13a) ), but also due to the prevalence of the [OlII]4959 line (see Section [A^ . The deviations from a Schechter function are likely explained by the mix 
of different lines contributing to the final LF, as the observed luminosity function is likely a combination of three different functions (for the three different 
emission lines). For higher luminosities, our data is in good agreement with the results from |Colbert et al"]^2Q13) and our slightly lower number densities 
are likely a result of accounting from the contribution of both [OlII]4959 and H/3 for the final volue probed. Most importantly, we compare our results with 
[Khostovan et al.H2Q15^ , who present the [OlII]+H/3 LF from HiZELS in COSMOS+UDS, and find a very good agreement. Right: [Oil] luminosity function 
at 2 = 2.2 from our SA22 survey and comparison with other lower redshift surveys such as|Gallego et al.|j2002|,|Takahashi et al.|^2007| andjBayliss et al.| 
j2011) . Similarly to the evolution of the Ho and [OlII]+H/3 luminosity functions, the [Oil] luminosity function shows a strong L* evolution. Such increase in 
L* with increasing redshift is even stronger by a few times than what is seen for other lines. On the other hand, there is little to no evidence of evolution in (p*. 
We also compare our results from SA22 with those from UDS and COSMOS (HiZELS) presented by |Khostovan et ar]j2015) and find a very good agreement, 
at both the bright and faint ends. 


For [On], in order to correct for dust extinction, we use the 
results from [Flayashi et al.| ( |2013^ that present a large study of 
[Oil] selected emitters at 2 ~ 1.47 and use Fla measurements for 
all of them (including stacking in Ha for the faintest [Oil] emit¬ 
ters) to calibrate [Oil], [Hayashi et al~] ( |2013^ shows that [Oil] se¬ 
lected samples at much lower luminosities than ours have typical 
dust extinctions of Ahc ~ 0.4, that their median extinction drops 
with increasing [Oil] luminosity, and that assuming Ana = 1.0 
for such emitters results in a significant overestimation of the real 
[Oil] luminosity. As a simple test of this conclusion (and its va¬ 
lidity at 2 = 2.2), we start by assuming Ahq = 1-0 for our 
sample of [Oil] emitters at 2 = 2.2, and find that we obtain 
PSFR = 1.0 ± 0.3M 0 yr~^ Mpc“^, about 4-5 times higher than 
that given by Ha ( Sobral^t^ y20]3a) . This suggests that the results 
presented by |Hayashi et al. j2013^ are also valid for 2 = 2.2. In¬ 
deed, if we use their results and instead extinction-correct the [Oil] 
luminosities assuming Ahc« ~ 0.2 (appropriate for our luminosity; 
this corresponds to Ajon] ~ 0.4), we obtain a good agreement with 
the value from Ha at the same redshift. However, we note that the 
uncertain faint-end slope of the [Oil] luminosity function at high 
redshift also plays a role. At lower redshift, it is found to be rela¬ 
tively shallow (a approximately in the range from —0.8 to —1.2), 
and thus here we fix it to a = —0.9, as this is the value found by 
the best studies at 2 ~ 1.5 ( |Ly et al.|2007[[Sobral et al.|2012^ . 


While e.g. |Ly et al.lpOOV) explored a double-blind [OlIl]-Ha 
survey (and found typical line ratios around ~ 1), no systematic 
investigation was performed in order to investigate the typical dust 
properties of [OlIl]-l-H/3 emitters at 2 > 1. However, as we have 


discussed, [Oll]-selected emitters are significantly less dusty than 
Ha selected line emitters (as an [Oll]-selection preferentially picks 
up dust-free sources and is biased against dusty sources). This is 
likely to be the case for [OlIl]-l-H/3 emitters, although we expect 
them to present typical dust extinction properties which are more 
similar to Ha-selected line emitters. However, given that we are us¬ 
ing a SFR calibration based directly on observed line fluxes which 
used the Ha line jLy et al.|2007| >, we will use the correction appro¬ 
priate for Ha-selected emitters. Ana = 1, corresponding to about 
1.5 mag at [OlIl]-bH/3. We obtain psFR = 0.14 ± O.O3M0yr“^ 
Mpc"^ 


Our results from Ha, [Olll] and [Oil] can be found in both 
Table and in Figure [T^ We find a very good agreement with 
the overall star formation history of the Universe (e.g.|Karim et al.| 
|201 l[[Sobral et al.|2013a| >, indicating that our simple assumptions 
for dust and AGN correction (or the combination of both) provide 
decent average estimates for the full population of galaxies (al¬ 
though it obviously fails on a source by source basis, and it will 
have strong dependences with mass, luminosity/etc - see e.g. |Ibar| 
let al.|2013|[Sobral et al.|20T4) . 

Our results confirm the rise of the star formation history up to 
2 = 2.2, but also indicate how important it is to apply appropriate 
dust extinction corrections jHayashi et al.pOO) that take into ac¬ 
count that [On] emitters are typically almost dust free. While there 
are significant uncertainties in the use of [Oil] as a star formation 
indicator for higher redshifts where Ha is not available, if a better 
understanding of the typical dust properties is obtained, and partic- 
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Redshift (z) 

Figure 13. Star formation history of the Universe using our CF-HiZELS 
SA22 data: 2 = 0.81 (Ha), .s = 1.4 ([Olll]+H/3) and .2 = 2.2 ([Oil]). Our 
results agree very well with the Ha star-formation history of the Universe 
(dashed line is the parameterisation from|Sobral et al.|2Q13a^. 


ularly if the faint-end slope is constrained, it should provide at least 
a competitive view when compared to e.g. the UV. 



Mpc (z = 0 . 81 ) 


Figure 14. Density of Ho emitters at 2 = 0.81 in our SA22 survey. The an¬ 
gular scale is converted to a physical scale. The smoothed distribution high¬ 
lights the large l-9cr over-density, which is approximately 15 Mpc across. 
It also shows the filamentary structure of the large scale structure in the 
Universe. 


5.4 Over-densities 

We explore our very wide area coverage over different fields to 
look for significant over-densities. In order to do this, we take two 
approaches: we smooth the distribution of sources, but also com¬ 
pute the 10-th nearest neighbour densities. As can be seen in Fig- 
ure[^ there appears to be a significant large-scale over-density of 
Ha emitters which contains ~ 300 candidate 2 = 0.81 Ha line 
emitters within a ~20-arcmin field (Figure [g. This includes a 
region where the number density of Ha emitters is ~ 10 times 
higher than the general field. This strong over-density of Ha emit¬ 
ters in the SA22 field has now been confirmed with the new KMOS 
instrument ( [Sobral et al.||2013b) . This ~ 8 ct over-density of Ha 
emitters is found within a volume of 3000 Mpc^ (co-moving). 

We note that this finding could be important to also interpret 
any over-density of star-forming galaxies detected at this or at even 
higher redshift. Without a careful analysis of the full galaxy pop¬ 
ulation, one could be misled to conclude that such structure could 
be a cluster and/or a proto-cluster. However, following|Sobral et al.| 
( |201 l^ , [Sobral et ai:] ( |2013b^ find that the structure presents galaxy 
number densities which are more typical of groups, or intermediate 
densities, but not clusters. In practice, the structure is likely a dense 
cosmic web structure or filamentary/wall-like as a group. |Darvish| 
|et al.| ( |2014^ find that the fraction of star-forming galaxies (Ha 
emitters) is much higher in filaments than in the general field or 
in clusters, making these easily detectable with wide Ha surveys. 

We also search for over-densities within the samples of [Olll] 
and [Oil] emitters, but find only mild over-densities (< 3 o). Nev¬ 
ertheless, these density fluctuations still result in significant cosmic 
variance, if only different parts of the fields would be investigated. 


5.5 Cosmic Variance and Cosmic convergence 

We explore the large samples spread over large, multiple areas 
(SA22, COSMOS, UDS) to empirically quantify cosmic variance 
affecting narrow-band surveys. In order to do this, we follow a 


similar procedure to |Sobral et al.|pbll| l, and divide the full sam¬ 
ple in several areas corresponding to individual pointing/cameras, 
roughly corresponding to 0.02 to 0.05 deg^. We then repeat the de¬ 
termination of the luminosity functions by randomly selecting these 
areas. We start by computing them for the smallest, contiguous 0.02 
deg^ areas and then go up by steps of 0.02, sampling these in ran¬ 
dom combinations within SA22, COSMOS and UDS. We first in¬ 
vestigate how a, 0* and L* are affected by sample variance, but 
find that, due to the depth of our data, a is in general completely 
unconstrained. Thus, for the following analysis, we fix a and con¬ 
centrate on studying (ft* and L* only. For each emission line, we 
obtain a total of a million realisations of the luminosity function, 
spread over the minimum to the maximum total area. Quoted errors 
(in %) are computed as the ratio between the standard deviation of 
a parameter for each bin in area (from all LF realisations with a 
given total area) and the actual parameter value. 

We note that our results can have a relatively general applica¬ 
tion for data-sets with similar depths, but that for data-sets that are 
e.g. significantly shallower, the errors measured with our method 
would become much larger, as they would be dominated by e.g. 
Poissonion errors and, in such cases (e.g. if depth is shallower than 
L*), (f* and L* become unconstrained, and thus the results of our 
study no longer apply (in order to lower the errors in the param¬ 
eters, deeper data is more efficient than probing larger areas). We 
also note that for [OlIl]-l-H/3 and [Oil] LFs, due to the depth of our 
data, we only probe down to ~ L* and thus, for the lowest vol¬ 
umes, our errors become completely Poissonion dominated, result¬ 
ing in very large variance which is due to low number statistics, not 
cosmic variance. Thus, we do not show the results for areas which 
are so small that low number statistics completely dominate. 

The results are presented in Figure [T^ As expected, we find 
that cosmic/sample variance is very strong for small volumes (< 
10® Mpc®), and strongly reduced with increasing area/volume. For 
the Ha luminosity function, for example, a total area of « 1 deg® 
(« 10® Mpc®) still results in significant scatter/significant errors 
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Figure 15. Top: The error/variance on (j>* , L* and SFRD (psfr) as func¬ 
tion of total probed co-moving volume, based on a million realisations, each 
sampling a random combination of UDS, COSMOS and SA22 for Ha LF 
at 2 ~ 0.8. The results show that the variance on all parameters/quantities 
decreases with total co-moving volume probed, but that it is significant for 
most volumes usually used in the literature. Our results are therefore ex¬ 
pected to have an error < 10%. Middle: The error/variance on (/>*, L* and 
SFRD (psfr) as function of total probed co-moving volume, based on a 
million realisations, each sampling a random combination of UDS, COS¬ 
MOS and SA22 for [Oiii]-hH/3 LF at 2 ~ 1.4. Bottom: The error/variance 
on <p*, L* and SFRD (psfr) us function of total probed co-moving vol¬ 
ume, based on a million realisations, each sampling a random combination 
of UDS, COSMOS and SA22 for [Oil] LF at 2 ~ 2.2. 


due to sample/cosmic variance, where L* is clearly the most af¬ 
fected parameter (up to >100% errors), followed by (j>* (~ 40 — 
50% errors). 

With our total area of ~ 9 deg^ the errors on (j)* are reduced 
to ~ 10 — 15%, while the errors on L* drop to ~ 8% and the error 
on SFRD drops to close to 5%. We therefore conclude that only by 
probing volumes comparable or larger than 10® Mpc® will the sam¬ 
ple variance fully drop below 10%. Our results are also very useful 
to interpret results from similar previous, current and future sur¬ 
veys of star-forming galaxies, particularly to interpret differences 
in number counts and luminosity functions, without having to rely 
on simulations and on other indirect methods/arguments. 

We find similar results for both [Oil] and [OlIl]-l-F[/3. We note 
that we use our results to obtain a better estimate of our errors 
by adding (in quadrature) the expected errors due to cosmic vari¬ 
ance, based on the total volume probed by each luminosity func¬ 
tion/measurement. This is a small contribution to the total error 
budget given that our survey probes a very large volume, but would 
be a significant amount to surveys probing small volumes. 


6 CONCLUSIONS 

We presented results from the largest contiguous narrow-band sur¬ 
vey in the near-infrared (J band). We surveyed «10deg^ of con¬ 
tiguous extragalactic sky and found a total of 5976 robust candidate 
emission-line galaxies (E > 3). By using deep CFHTLS ugriz 
and UKIDSS DXS J and K, a large sample of available spectro¬ 
scopic redshifts from VVDS and VIPERS from the literature, and 
by obtaining new spectroscopic follow-up observations with MOS- 
FIRE and EMOS, and by combining our sample with HiZELS, we 
derive the largest samples of emission-line selected galaxies. We 
find by far the largest sample of Ha emitters at 2 ~ 0.8: 3471 
sources (~ 400 spectroscopically confirmed), obtaining the most 
accurate measurement of the star formation rate density at that cos¬ 
mic epoch. We also present the largest sample (1341) of [OlIl]/H/3 
emitters at 2 ~ 1.4 and present the first large sample (572 sources) 
of [Oil] emitters at the peak of the star formation history (2 = 2.2). 
Our main conclusions are: 

• Our large spectroscopic sample from FMOS and MOSFIRE 
allows us to confirm that the [Nll]/Ha correction as a func¬ 
tion of EW(Ha-l-[Nll]) from SDSS is applicable to 2 = 0.81 
with no evolution. We nonetheless provide a simpler linear fit, 
based on our data, that can be used to correct similar data-sets: 
/([NIIJ/Ha) = —0.296xlogio(EWHc«+[Nii])+0.8. Ourspectro- 
scopic Ha sample also shows that our Ha emitters have a metallic- 
ity of 12 + logic (0/H) = 8.56 ±0.05, slightly sub-solar. We also 
make significant detections of the [Sll] doublet, with a line ratio 
of I([SII]6716)/I([SII]6731) = 1.33 ± 0.08, implying an electron 
density of 40-200cm"®. The [Sll]6716/I(Ha)= 0.14 ± 0.02 ratio 
also implies ionisation potential of logioU = —3.9 ± 0.5 cm“^. 

• We obtain Ha, [OlIl]-l-H/? and [Oil] luminosity functions at 
2 = 0.8,1.4, 2.2 with the largest statistical samples, reaching up to 
the highest luminosities. We find a strong luminosity evolution (L*) 
in the luminosity function of all the lines with increasing redshift up 
to at least 2 ~ 2.2, with a less significant (but present) (j>* evolution. 
This is consistent with the evolution seen across redshift presented 
by |Khostovan et al!] ( |2015| l. 

• We show that the [OlIl]-l-H(? luminosity function at 2 = 1.4 is 
very hard to interpret in general, particularly due to the complicated 
contribution from the two different [Olll] lines, and, to a lesser ex¬ 
tent, to the contribution of H/3. While we find that [Olll]5007 (only) 
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Table 5. The luminosity function and star formation rate density evolution for 2 : = 0.8,1.4, 2.2. L*, 0* and a. values are derived without any correction for 
dust extinction. The measurements connected for 1 mag extinction at Hq for Ha and [OlII]+H/3 emitters (as our SFR calibration is based on observed Ha 
and [OlII] line fluxes, thus the correction for Ha emitters is the most appropriate), but follow [Hay ashi et al.H2013} for [Oil] selected emitters (Ana = 0.2). 
Columns present the redshift, break of the luminosity function, normalisation faint-end slope (a) of the luminosity functions. The two right 

columns present the star formation rate density at each redshift based on integrating the luminosity function down to each observational limit (41.4 for Ha, 
42.1 for [OlII] and 42.6 for [Oil] in logL) and for a full integration, and include dust connections. The two columns immediately to the left present similar 
measurements but for luminosity densities. Star formation rate densities include a correction for AGN contamination of 10% at 2 : = 0.8 for Ha (see jGarnj 
jet al.|20lQ) , up to L*, and a correction given by 0.59 x logjQ(L/L*) + 0.112 for L > L* jSobral et al.|2015') . For [OlII]+H^ at 2 : ~ 1.4, we assume an 
AGN contamination of 15%, while we assume a contamination of 20% at z = 2.2 for [Oil] emitters. We only fit a for the Ha LF and fix it for all the other 
luminosity functions. L* and <j>* are obtained by fixing a = —1.6 (Ha), and the 1 cr errors on L* and <j>* are derived from such fits (with fixed a). 
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dominates the sample of spectroscopically confirmed [Olll]+H/3 
emitters, representing ~ 50% of the full sample, [Olll]4959 (only) 
represents 27%, almost twice as common as H/3. We find that 7% 
of the [Olll]+H/3 sources are actually detections of both [Olll] lines 
at opposite wings of the filter. 

• We present the first [Oil] luminosity function at 2 = 2.2, and 
find a very strong evolution from 2 ~ 0 to 2 ~ 2.2, much stronger 
than that seen for either Ha or [Olll]+H/3. By correcting for dust 
extinction using|Hayashi et al.|(|2013|l the star formation rate den¬ 
sity based on [Oil] is in excellent agreement with Hq jSobral et al.| 
|2013a^ . If 1 magnitude of extinction at Ha was used for [Oil] emit¬ 

ters instead, the star formation rate density would have been over¬ 
estimated by a factor of 2-3. Thus, if the |Hayashi et al.|f2013^ cal¬ 
ibration is used (and if it remains valid for even higher redshift), 
[Oil] may be a reasonably good way to measure star formation rate 
density beyond 2 ~ 2.2, as Ha gets redshifted out of the K band. 

• We find a reasonable good agreement between the star forma¬ 
tion rate density from [Olll] and that from Ha at a similar redshift. 
We nonetheless caution that without a detailed investigation into 
the nature of [Olll]-l-H/l emitters (and without robustly separating 
them from H/3 emitters, interesting on their own), using [Olll]-l-H/3 
as a star formation indicator at high redshift is highly unreliable. 

• For the planned Euclid wide survey ( [Laureijs et al.|20i^ flux 
limit of > 3 X 10“^®ergs“^, our sample is dominated by Ha 
emitters (67%), followed by [Olll]-bH/3 emitters (12%), but still 
with some [Oil] emitters (3%); the remaining 18% are rarer emit¬ 
ters spread over a wide range of redshifts. The raw (observed) 
density of Ha emitters with > 3 X 10 '^®ergs ^ at 2 ~ 0.8 
(A ~ 1.2 fim) is 8.66xl0~® Mpc“®, while at the same wave¬ 
length the observed number density of [OlIl]-bH/3 emitters is about 
5 times lower (1.6x 10“® Mpc~®), and the number density of [Oil] 
emitters (0.3x 10“® Mpc~®) is almost 30 times lower than Hq at 
A ~ 1.2 fim. Our fully spectroscopically confirmed sample con¬ 
firms these numbers, predicting that ~ 70 — 80% of all line emitters 
found at A ~ 1.2 fim will be Ha. 

• We find significant over-densities, out of which the strongest 
one is found at 2 = 0.8 and traced by Ha emitters. It is an 8.5 a 
(confirmed with KMOS) over-density of Hq emitters, where the 
number density is a factor ~ 10 x higher than the average; this 
is consistent with group-like densities, and most likely a rich fila¬ 
mentary structure, similar to what has been studied in|Darvish et al.| 
(|2014). Only mild over-densities are found for [Olll]-l-H/l and [Oil] 


emitters, although such over-densities would be much harder to find 
given the relatively shallow data when compared to Hq. 

• We take advantage of the large volumes/area and multi¬ 
ple fields to sub-divide the samples in randomised areas and 
provide a robust empirical measurement of sample/cosmic vari¬ 
ance for the different lines/redshifts. We find that surveys for 
star-forming/emission-line galaxies can only overcome cosmic- 
variance (errors < 10%) if they are based on volumes > 5 x 
10®Mpc®. In other words, multiple/different volumes adding up 
to < 5 X 10® Mpc® show variance which results in errors being 
> 10%. 

• Errors due to sample (cosmic) variance on surveys probing 
~ lO'^ Mpc® and ~ 10® Mpc® are typically very high: ~ 300% 
and ~ 40 — 60%, respectively. Focusing only on L* and (j)*, the 
latter is the most affected parameter for large volumes, while L* is 
completely un-determined for volumes < 10® Mpc®. 
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APPENDIX A: CATALOGUE OE CE-HIZELS 
NARROW-BAND EMITTERS 

The catalogue of narrow-band emitters in SA22 is presented in Ta¬ 
ble A.l. It contains IDs, Right Ascension (RA, J2000), Declination 
(Dec, J2000), narrow-band magnitude (NB, AB), broad-band mag¬ 
nitude (J, AB), the significance of the narrow-band excess (S, esti¬ 
mated in 2" apertures), estimated flux (logj^p), estimated observed 
EW (A), and a flag for those that are classified as Ha at a = 0.8 
(1), [OIII] or H/3 at z ~ 1.4 (2) and [Oil] at ^ = 2.20 (3). Unclas¬ 
sified sources as flagged with 0 and candidate stars are identified 
with flag -1. We note that there are a few (15) sources with some 
extreme EWs (observed EWs > 10'* A [observed]) and 17 sources 


Table Bl. Hq Luminosity Function at z = 0.81 from our SA22 survey. 
Note that Lhci has been corrected for [Nil] contamination, but not for dust 
extinction. The total volume per bin is 8.6 X 10® Mpc®, assuming a top hat 
filter, but 0 corr has been corrected for both incompleteness and the fact 
that the filter profile is not a perfect top hat. 


logLHa 

Sources 

0 obs 

0 corr 

z = 

: 0.81 


# 


Mpc 

3 

Mpc 

3 

41.40 

± 

0.025 

445 

± 

21 

-1.98 

± 

0.02 

-2.01 

± 

0.02 

41.45 

± 

0.025 

433 

± 

20 

-1.99 

± 

0.02 

-2.03 

± 

0.02 

41.50 

± 

0.025 

370 

± 

19 

-2.06 

± 

0.02 

-2.09 

± 

0.02 

41.55 

± 

0.025 

322 

± 

17 

-2.12 

± 

0.02 

-2.15 

± 

0.03 

41.60 

± 

0.025 

268 

± 

16 

-2.20 

± 

0.03 

-2.23 

± 

0.03 

41.65 

± 

0.025 

248 

± 

15 

-2.24 

± 

0.03 

-2.27 

± 

0.03 

41.70 

± 

0.025 

201 

± 

14 

-2.33 

± 

0.03 

-2.35 

± 

0.03 

41.75 

± 

0.025 

159 

± 

12 

-2.43 

± 

0.04 

-2.46 

± 

0.04 

41.80 

± 

0.025 

109 

± 

10 

-2.59 

± 

0.04 

-2.63 

± 

0.05 

41.85 

± 

0.025 

76 

± 

8 

-2.75 

± 

0.05 

-2.79 

± 

0.06 

41.90 

± 

0.025 

73 

± 

8 

-2.77 

± 

0.05 

-2.81 

± 

0.06 

41.95 

± 

0.025 

47 

± 

6 

-2.95 

± 

0.07 

-2.98 

± 

0.07 

42.00 

± 

0.025 

46 

± 

6 

-2.96 

± 

0.07 

-3.01 

± 

0.07 

42.05 

± 

0.025 

23 

± 

4 

-3.26 

± 

0.10 

-3.31 

± 

0.10 

42.10 

± 

0.025 

21 

± 

4 

-3.29 

± 

0.10 

-3.33 

± 

0.11 

42.15 

± 

0.025 

13 

± 

3 

-3.50 

± 

0.14 

-3.54 

± 

0.14 

42.20 

± 

0.025 

9: 

±: 

3 

-3.64 

± 

0.17 

-3.68 

± 

0.17 

42.25 

± 

0.025 

5: 

±: 

2 

-3.93 

± 

0.26 

-4.00 

± 

0.26 

42.30 

± 

0.025 

5: 

±: 

2 

-3.93 

± 

0.26 

-3.93 

± 

0.26 

42.45 

± 

0.025 

13 

± 

3 

-4.44 

± 

0.14 

-4.44 

± 

0.14 

42.65 

± 

0.025 

3: 

± 

1 

-5.08 

± 

0.33 

-5.02 

± 

0.33 


Table B2. [OIII]+H/3 Luminosity Function - this has not been corrected for 
extinction. The total volume per bin is 24.8 x 10^ Mpc^, assuming a top 
hat filter, but (p corr has been corrected for both incompleteness and the fact 
that the filter profile is not a perfect top hat. 

logL[oiii] 

Sources 

0 obs 

0 eorr 

z = 1.4 

# 

Mpc“® 

Mpc“® 

42.07 ±0.03 

285 ± 16 

-2.72 ± 0.03 

-2.71 ± 0.04 

42.13 ±0.03 

283 ± 16 

-2.72 ± 0.03 

-2.71 ± 0.04 

42.19 ±0.03 

187 ± 13 

-2.90 ±0.03 

-2.91 ± 0.04 

42.25 ± 0.03 

104 ± 10 

-3.16 ±0.04 

-3.12 ±0.06 

42.31 ± 0.03 

74 ± 8 

-3.30 ±0.05 

-3.21 ± 0.08 

42.37 ±0.03 

32 ±5 

-3.67 ± 0.08 

-3.53 ±0.15 

42.43 ±0.03 

23 ±4 

-3.81 ± 0.10 

-3.63 ±0.19 

42.49 ±0.03 

17 ±4 

-3.94 ±0.12 

-3.78 ±0.23 

42.55 ± 0.03 

8± 2 

-4.27 ± 0.19 

-4.10 ± 0.43 

42.61 ± 0.03 

6± 2 

-4.39 ±0.23 

-4.21 ± 0.56 

42.73 ±0.03 

14 ± 3 

-4.03 ±0.14 

-4.28 ± 0.29 

42.85 ± 0.09 

4± 2 

-4.57 ± 0.30 

-5.03 ±0.53 

43.03 ±0.09 

1± 1 

-5.17 ± 0.38 

-5.63 ±0.38 


with E > 100, but that these may be supernovae and/or strongly 
variable sources (see e.g. |Matthee et al.|2015) , apart from real ex¬ 
tremely rare sources. Note that only the online version contains the 
full catalogue - here only three entries of the table are shown as 
examples of the entire catalogue. 
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Table Al. Example entries form the catalogue of all S > 3 narrow-band sources selected in the SA22 field from this paper. The full catalogue is available 
on-line. 


ID 

R.A. 

(J2000) 

Dec. 

(12000) 

NBj 

(AB) 

J 

(AB) 

S 

log Flux 
ergs“^ cm~^ 

EWob, 

A 

Class. Flag 

CFHIZELS-SM14-99905 

22 09 11.64 

+01 23 13.25 

20.43+0.06 

21.40+0.06 

10.1 

-15.48 

170.8 

1 

CFHIZELS-SM14-98588 

22 19 34.55 

+00 24 51.36 

20.41+0.06 

22.24+0.11 

13.2 

-15.33 

657.2 

2 

CFHIZELS-SM14-97547 

22 11 31.53 

-00 37 16.69 

20.40+0.07 

20.82+0.05 

5.1 

-15.72 

52.4 

3 


Table B3. [On] Luminosity Function not corrected for dust extinction. The 
total volume per bin is 26.2 x 10^ Mpc^, assuming a top hat filter, but (p 
corr has been corrected for both incompleteness and the fact that the filter 
profile is not a perfect top hat. 


logLfoii] 

Sources 

</> obs 

(p corr 

z = 

= 2.2 

# 

Mpc 

3 

Mpc 

3 

42.60 

+ 

0.025 

141 + 11 

-2.97 

+ 

0.04 

-2.88 

+ 

0.11 

42.65 

+ 

0.025 

118+ 10 

-3.04 

+ 

0.04 

-2.96 

+ 

0.06 

42.70 

+ 

0.025 

73 + 8 

-3.25 

+ 

0.05 

-3.14 

+ 

0.12 

42.75 

+ 

0.025 

53 + 7 

-3.39 

+ 

0.06 

-3.31 

+ 

0.11 

42.80 

+ 

0.025 

30 + 5 

-3.64 

+ 

0.09 

-3.54 

+ 

0.14 

42.85 

+ 

0.025 

18 + 4 

-3.86 

+ 

0.12 

-3.75 

+ 

0.15 

42.90 

+ 

0.025 

15 + 3 

-3.94 

+ 

0.13 

-3.84 

+ 

0.16 

42.95 

+ 

0.025 

6 + 2 

-4.34 

+ 

0.23 

-4.21 

+ 

0.25 

43.00 

+ 

0.025 

2+1 

-4.82 

+ 

0.53 

-4.69 

+ 

0.58 

43.10 

+ 

0.075 

5 + 2 

-4.99 

+ 

0.26 

-4.99 

+ 

0.26 

43.25 

+ 

0.075 

2+1 

-5.52 

+ 

0.53 

-5.52 

+ 

0.53 


APPENDIX B: LUMINOSITY FUNCTIONS 
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